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Abstract 

We give a pedagogical introduction to topologically ordered states of matter, with the aim of 
familiarizing the reader with their axiomatic topological quantum field theory description. We in¬ 
troduce basic noninteracting topological phases of matter protected by symmetries, including the 
Su-Schrieffer-Heeger model and the one-dimensional p -wave superconductor. The defining proper¬ 
ties of topologically ordered states are illustrated explicitly using the toric code and - on a more 
abstract level - Kitaev’s 16-fold classification of two-dimensional topological superconductors. Sub¬ 
sequently, we present a short review of category theory as an axiomatic description of topological 
order in two-dimensions. Equipped with this structure, we revisit Kitaev’s 16-fold way. 
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cations”, 25-29 May 2015, Universidad de los Andes, Bogota, Colombia 


1 



CONTENTS 


I. Introduction to topological phases in condensed matter 4 

A. The notion of topology 4 

B. Classification of noninteracting fermion Hamiltonians: The 10-fold way 7 

1. Classification with respect to time-reversal and particle-hole symmetry 7 

2. Flatband Hamiltonians and homotopy groups 9 

3. Topological invariants 12 

C. The Su-Schrieffer-Heeger model 16 

D. The one-dimensional p- wave superconductor 18 

E. Reduction of the 10-fold way classification by interactions: Z — > Z 8 in class BDI 22 

II. Examples of topological order 26 

A. The toric code 26 

1. Ground states 28 

2. Topological excitations 31 

B. The two-dimensional p-wave superconductor 34 

1. Argument for the existence of Majorana bound states on vortices 40 

2. Bound states on vortices in two-dimensional chiral p- wave superconductors 42 

3. Non-Abelian statistics of vortices in chiral p-wave superconductors 43 

4. The 16-fold way 46 

III. Category theory 49 

A. Fusion Category 49 

1. Diagrammatics 51 

2. F-moves and the pentagon equation 52 

3. Gauge freedom and its fixing 54 

4. Quantum dimensions and Frobenius Schur indicators 55 

5. Examples 57 

B. Braiding Category 58 

1. Topological spin 59 

2. Ribbon equation 61 

3. Vafa’s Theorem 62 


2 


C. Modular matrices 62 

1. The S matrix 63 

2. Verlinde Formula 64 

3. Obstruction for theories with multiplicities 66 

4. The T matrix 67 

D. Examples: The 16-fold way revisited 67 

1. Case: C (1 ) odd 67 

2. Case C^ 1 ) = 2 mod 4 68 

3. Case C^ 1 ) = 0 mod 4 69 

Acknowledgements 70 

References 70 


3 


I. INTRODUCTION TO TOPOLOGICAL PHASES IN CONDENSED MATTER 


A. The notion of topology 

In these lectures we will learn how to categorize and characterize some phases of matter 
that have topological attributes. A topological property of a phase, such as boundary modes 
(in an open geometry), topological response functions, or the character of its excitations, is 
described by a set of quantized numbers, related to so-called topological invariants of the 
phase. The quantization immediately implies that topological properties are universal (they 
can be used to label the topological phase) and in some sense protected, because they cannot 
change smoothly when infinitesimal perturbations are added. Topological properties, in the 
sense that we want to discuss them here, can only be defined for 

• spectrally gapped ground states on a manifold without boundary of 

• local Hamiltonians at 

• zero temperature. 

The spectral gap allows to define an equivalence class of states, i.e., a phase, with the help 
of the adiabatic theorem. Two gapped ground states are in the same phase if there exists 
an adiabatic interpolation between their respective Hamiltonians, such that the spectral gap 
above the ground state as well as the locality is preserved for all Hamiltonians along the 
interpolation. 

Often it is useful to further modify these rules to define topological phases that are 
subject to symmetry constraints. We refer to topological states as being protected/enriched 
by a symmetry group G, if the Hamiltonian has a symmetry G and only G-preserving 
interpolations are allowed. Since the G-preserving interpolations are a subset of all local 
interpolations, it is clear that symmetries make a topological classification of Hamiltonians 
more refined. 

The locality of a Hamiltonian is required to guarantee the quantization of topological 
response functions and to distinguish topological characterizations depending on the di¬ 
mensionality of space. If we were not to impose locality, any system could in essence be 
zero-dimensional and there would be no notion of boundary states (which arc localized over 
short distances) or point-like and line-like excitations etc. 
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Equipped with this definition of a topological phase, the exploration of topological states 
of matter above all poses a classification problem. We would like to know how many phases 
of quantum systems exists, that can be distinguished by their topological properties. We 
would like to obtain such a classification while imposing any symmetry G that is physically 
relevant, such as time-reversal symmetry, space-group or point-group symmetries of a crystal, 
particle-number conservation etc. To identify the right mathematical tools that allow for 
such a classification and to guarantee its completeness is a subject of ongoing research. 
Here, we shall focus on aspects of this classification problem, which are well established and 
understood. 

Most fundamental is a distinction between two types of topological states of matter: 
Those with intrinsic (long-range entangled) topological order 1 and those without. This 
notion is also core to the structure of these lecture notes. In this Section, we only discuss 
phases without intrinsic topological order, while the ensuing two Sections are devoted to 
states with intrinsic topological order. A definition of intrinsic topological order can be 
based on several equivalent characterizations of such a phase, of which we give three: 

• Topological ground state degeneracy: On a manifold without boundary, the degeneracy 
of gapped topologically degenerate ground states depends on the topological properties 
of the manifold. There are no topologically degenerate ground states if the system is 
defined on a sphere. The matrix elements of any local operator taken between two 
distinct topologically degenerate ground states vanishes. 

• Fractionalized excitations: There exist low-energy excitations which are point-like 
[in two dimensions (2D) or above] or line-like [in three dimensions (3D) or above]. 
These excitations carry a fractional quantum numbers as compared to the microscopic 
degrees of freedom that enter the Hamiltonian (for example, a fractional charge), are 
deconfined and dynamical (i.e, free to move in the low-energy excited states). 

• Topological entanglement entropy: The entanglement entropy between two parts of a 
system that is in a gapped zero-temperature ground state typically scales with the size 
of the line/surface that separates the two regions (“area-law entanglement”). Topo¬ 
logically ordered, long-range entangled states have a universal subleading correction 
to this scaling that is characteristic for the type of topological order. 
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(Note that these statements, as many universal properties we discuss, are only strictly 
true in the thermodynamic limit of infinite system size. For example, in a finite system, 
the ground state degeneracy is lifted by an amount that scales exponential in the system 
size.) As fractionalized excitations in the above sense may only exist in two or higher 
dimensions, intrinsic topological order cannot be found in one-dimensional (ID) phases of 
matter. Further, for intrinsic topological order to occur, interactions are needed in the 
system. 

Examples of topologically nontrivial phases (both with and without intrinsic topological 
order) exist in absence of any symmetry. However, most of the phases without intrinsic 
topological order belong to the so-called symmetry protected topological (SPT) phases. In 
these cases, the topology is protected by a symmetry. These phases almost always possess 
topologically protected boundary modes when defined on a manifold with boundary, except 
if the boundary itself breaks the protecting symmetry (as could be the case with inversion 
symmetry, for example). 

In contrast, phases with intrinsic topological order are not necessarily equipped with 
boundary modes, even if the boundary of the manifold preserves the defining symmetries of 
the phase. If the definition of a phase with intrinsic topological order relies on symmetries, 
it is named symmetry enriched topological phase (SET). 

An alternative characterization of topological properties of a phase uses the entanglement 
between different subsystems. While we opt not to touch upon this concept here, we want 
to make contact to the ensuing terminology: All phases with intrinsic topological order are 
called long-range entangled (LRE). The term short-range entangled (SRE) phase is often 
used synonymously with “no intrinsic topological order”. (Some authors count 2D phases 
with nonvanishing thermal Hall conductivity, such as the p + ip superconductors, but no 
intrinsic topological order unless gauged, as LRE.) 

In these lecture notes, we will encounter two classifications of a subset of topological 
phases. The following Subsection introduces the complete classification of non-interacting 
fermionic Hamiltonians with certain symmetries (which have no intrinsic topological order). 
Section III is concerned with the unified description of 2D phases with intrinsic topological 
order in absence of any symmetries. 
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B. Classification of noninteracting fermion Hamiltonians: The 10-fold way 

We have stated that SPT order in SRE states manifests itself via the presence of gapless 
boundary states in an open geometry. In fact, there exists a intimate connection between 
the topological character of the gapped bulk state and its boundary modes. The latter are 
protected against local perturbations on the boundary that (i) preserve the bulk symmetry 
and (ii) induce no intrinsic topological order or spontaneous symmetry breaking in the 
boundary modes. This bulk-boundary correspondence can be used to classify SPT phases. 
Two short-range entangled phases with the same symmetries belong to a different topological 
class, if the interface between the two phases hosts a state in the bulk gap and this state 
cannot be moved into the continuum of excited states by any local perturbation that obey 
(i) and (ii). Equivalently, to change the topological attribute of a gapped bulk state via any 
smooth changes in the Hamiltonian, the bulk energy gap has to close and reopen. 

In Ref. 2 Schnyder et al. use this bulk-boundary correspondence to classify all nonin¬ 
teracting fermionic Hamiltonians. For the topological phases that they discuss, two funda¬ 
mental symmetries, particle-hole symmetry (PHS) and time-reversal symmetry (TRS), are 
considered. In the following, we will review the essential results of this classification W 1 


1. Classification with respect to time-reversal and particle-hole symmetry 

Symmetries in quantum mechanics are operators that have to preserve the absolute value 
of the scalar product of any two vectors in the Hilbert space. They can thus be either unitary 
operators, preserving the scalar product, or antiunitary operators, turning the scalar product 
into its complex conjugate (up to a phase). For a unitary operator to be a symmetry of a 
given Hamiltonian H , the operator has to commute with H. Consequently, the Hamiltonian 
can be block diagonalized, where each block acts on one eigenspace of the unitary symmetry. 
If H has a unitary symmetry, we block-diagonalize it and then consider the topological 
properties of each block individually. This way, we do not have to include unitary symmetries 
(except for the product of TRS and PHS and the omnipresent particle number conservation) 
in the further considerations, as we will not focus on the burgeoning field of crystalline 
topological insulators. 

A fundamental antiunitary operator in quantum mechanics is the reversal of time T. Let 
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us begin by recalling its elementary properties. If a given Hamiltonian H is TRS, that is, 


THT- 1 = +H, 


(la) 


the time-evolution operator at time t should be mapped to the time-evolution operator at 
—t by the operator T 


T e- itH T~ l =e- TiT ~ ltH 
= e-^ H . 


(lb) 


We conclude that the reversal of time is indeed an antiunitary operator TiT^ 1 = — i. It 
can be represented as T = TIC, where /C denotes complex conjugation and T is a unitary 
operator. Applying the reversal of time twice on any state must return the same state up 
to an overall phase factor e 1 ^* 

e '= T 2 = T(T t ) -1 => T — e^T J , T J = t&T. (lc) 


Inserting the two last equations into one another, one obtains T = e 2l< ^T, i.e., e 21 ^ has to 
equal +1. We conclude that the time-reversal operator either squares to +1 or to —1 


T 2 = +1, T 2 = -1. 


(Id) 


The second fundamental antiunitary symmetry considered here is charge conjugation V. 
Its most important incarnation in solid state physics is found in the theory of supercon¬ 
ductivity. In an Andreev reflection process, an electron-like quasi particle that enters a 
superconductor is reflected as a hole-like quasi particle. The charge difference between inci¬ 
dent and reflected state is accounted for by adding one Cooper pair to the superconducting 
condensate. In the mean-held theory of superconductivity, the energies of the electron-like 
state and the hole-like state are equal in magnitude and have opposite sign, giving rise to 
the PHS. In this case, rather than being a fundamental physical symmetry of the system 
like TRS is, PHS emerges clue to a redundancy in the mean-held description. We define a 
(single-particle) Hamiltonian H to be PHS if 

VHV- 1 = +H. (le) 

In order to also reverse the sign of charge, V has to turn the minimal coupling p — i eA 
into p + i eA, where p is the momentum operator and A is the electromagnetic gauge 



potential. This is achieved by demanding ViV = —i. We conclude that V is indeed an 
antiunitary operator that can be decomposed as V = PIC, where P is a unitary operator. As 
a consequence, the reasoning of Eq. (lc) also applies to V and we conclude that the charge 
conjugation operator either squares to +1 or to —1 

V 2 = +1, V 2 = -1. (If) 

In the case where the operators T and V are both symmetries of H, their product is 
also a symmetry of H. We call this product chiral transformation C := TV. It is a unitary 
operator. The Hamiltonian p[ transforms under the chiral symmetry as 

CHC~ l = +H. (lg) 

(It is important to note that both V and C anticommute rather than commute with the 
single-particle first-quantized Hamiltonian 'H a ,a 1 that we will introduce below.) Observe that 
a Hamiltonian can have a chiral symmetry, even if it possesses neither of PHS and TRS. We 
can now enumerate all combinations of the symmetries V, T, and C that a Hamiltonian can 
obey, accounting for the different signs of T 2 and V 2 . There are in total ten such symmetry 
classes, listed in Tab. I. The main result of Schnyder et al. in Ref. 2 is to establish how 
many distinct phases with protected edge modes exist on the (d — l)-dimensional boundary 
of a phase in d dimensions. We find three possible cases: If there is only one (topologically 
trivial) phase, the entry 0 is found in Tab. I. If there are exactly two distinct phases (one 
trivial and one topological phase), Z 2 is listed. Finally, if there exists a distinct topological 
phase for every integer, Z is listed. 

2. Flatband Hamiltonians and homotopy groups 

There are several approaches to obtain the entries Z 2 and Z in Tab. I. For one, the 
theory of Anderson localization can be employed to determine in which spatial dimensions 
boundaries can host localization-protected states (the topological surface states) under a 
given symmetry. This was done by Schnyder et al. in Ref. 2. Kitaev, on the other hand, 
derived the table using the algebraic structure of Clifford algebras in the various dimensions 
and symmetry classes. 1 In mathematics, this goes under the name K-theory. 

Here, we want to give a flavor of the mathematical structure behind the table by con¬ 
sidering two examples. To keep matters simple, we shall restrict ourselves to the situation 
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TABLE I. Symmetry classes of noninteracting fermionic Hamiltonians from Refs. 3 and 4. The 
columns contain from left to right: Cartan’s name for the symmetry class; the square of the time 
reversal operator, the particle-hole operator, and the chiral operator (0 means the symmetry is 
not present); the group of topological phases that a Hamiltonian with the respective symmetry 
can belong to for the dimensions d = 1, • • • ,8 of space. The first two rows are called “complex 
classes”, while the lower eight rows are the “real classes”. The homotopy groups of the former show 
a periodicity with period 2 in d, while those of the latter have a period 8 in d (Bott periodicity). 
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where the system is translationally invariant and periodic boundary conditions are imposed. 
In second quantization, the Hamiltonian H has the Bloch representation 

H= I d d k^l(k)n a , a '{k)i> a ,(k), (2a) 

where ^ a {k) creates a fermion of flavor a = 1, • • • , N at momentum k in the Brillouin zone 
(BZ) and the summation over a and a! is implicit. The flavor index may represent orbital, 
spin, or sublattice degrees of freedom. Energy bands are obtained by diagonalizing the 
N x N matrix 'H(k) at every momentum k G BZ with the aid of a unitary transformation 
U(k) 

U^{k)U{k)U(k) = diag [£ m+n {k), ■ ■ • , e n+1 (k), e n {k), • • • ^(fc)] , ( 2b ) 

where the energies are arranged in descending order on the righthand side and n, m € Z 
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such that n + m = N. So as to start from an insulating ground state, we assume that there 
exists an energy gap between the bands n and n + 1 and that the chemical potential // lies 
in this gap 

£ n (k) < fi < e n+1 (k), Vfc G BZ. (3) 


The presence of the gap allows us to adiabatically deform the Bloch Hamiltonian 'H(k) to 
the hatband Hamiltonian 


Q(k) := U(k) 



(4a) 


that assigns the energy —1 and +1 to all states in the bands below and above the gap, respec¬ 
tively. This deformation preserves the eigenstates, but removes the nonuniversal information 
about energy bands from the Hamiltonian. 

In other words, the degenerate eigenspaces of the eigenvalues ±1 of Q(k) reflect the par¬ 
titioning of the single-particle Hilbert space introduced by the spectral gap in the spectrum 
of 'H(k). The degeneracy of its eigenspaces equips Q(k) with an extra U(n) x U(m) gauge 
symmetry: While the (n + m) x (n + m) matrix U(k) of Bloch eigenvectors that diagonalizes 
Q(k) is an element of U(n + m) for every k € BZ, we are free to change the basis for its 
lower and upper bands by a U(n) and U(m) transformation, respectively. Hence Q(k) is an 
element of the space C 0 := U(n + m)/[XJ(n) x U(m)] defining a map 


Q : BZ^C 0 . 


( 5 ) 


The group of topologically distinct maps Q, or, equivalently, the number of topologically 
distinct Hamiltonians +L, is given by the homotopy group 


Ti (Co) ( 6 ) 

for any dimension d of the BZ. (The homotopy group is the group of equivalence classes of 
maps from the d -dimensional sphere to a target space, in this case Co- Even though the BZ 
is a d-dimensional torus, it turns out that this difference between torus and sphere does not 
affect the classification as discussed here.) 

For example, in d = 2 we have n 2 (C 0 ) = Z. A physical example of a family of Hamil¬ 
tonians that exhausts the topological sectors of this group is found in the integer quantum 
Hall effect. The incompressible ground state with r £ N filled Landau levels is topologically 
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distinct from the ground state with N 3 r' ^ r filled Landau levels. Two different patches of 
space with r and r' filled Landau levels have |r — r'\ gapless edge modes running at their in¬ 
terface, reflecting the bulk-boundary correspondence of the topological phases. In contrast, 
7 r 3 (Cq) = renders all noninteracting fermionic Hamiltonians in 3D space topologically 
equivalent to the vacuum, if no further symmetries besides the U(l) charge conservation are 
imposed. 

As a second example, let us discuss a Hamiltonian that has only chiral symmetry and 
hence belongs to the symmetry class AIII. The chiral symmetry implies a spectral symmetry 
of 'H(k). If gapped, 'H(k) must have an even number of bands N = 2 n, n € Z. When 
represented in the eigenbasis of the chiral symmetry operator C, the spectrally flattened 
Hamiltonian Q(k) and the chiral symmetry operator have the representations 


Q(k) 


( 0 q(k)\ 

Vz f (fc) 0 ) 



(7a) 


respectively. From Q{k ) 2 = 1, one concludes that q{k) can be an arbitrary unitary matrix. 
We are thus led to consider the homotopy group it rf (C 1 ) of the mapping 


q: BZ^C 1= U(n). 


(7b) 


For example, in d, — 1 spatial dimensions ^(C^) = Z. A tight-binding model with non-trivial 
topology that belongs to this symmetry class will be discussed in Sec. IC. 

With these examples, we have discussed the two complex classes A and AIII. In the real 
classes, which have at least one antiunitary symmetry, it is harder to obtain the constraints 
on the spectrally flattened Hamiltonian Q(k). The origin for this complication is that the 
antiunitary operators representing time-reversal and particle-hole symmetry relates Q(k) 
and Q(—k ) rather than acting locally in momentum space. 


3. Topological invariants 

Given a gapped noninteracting fermionic Hamiltonian with certain symmetry properties 
in d-dimensional space, one can use Tab. I to conclude whether the system can potentially 
be in a topological phase. However, to understand in which topological sector the system is, 
we have to do more work. To obtain this information, one computes topological invariants 
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or topological quantum numbers of the ground state. Such invariants are automatically 
numbers in the group of possible topological phases (Z or Z 2 ). For many of them, a variety 
of different-looking but equivalent representations are known. 

To give concrete examples, we shall discuss the invariants for all Z topological phases 
found in Tab. I. These are called Chern numbers in the symmetry classes without chiral 
symmetry and winding numbers in the classes with chiral symmetry. 

In physics, topological attributes refer to global properties of a physical system that is 
made out of local degrees of freedom and might only have local, i.e., short-ranged, correla¬ 
tions. The distinction between global and local properties parallels the distinction between 
topology and geometry in mathematics, where the former refers to global structure, while 
the latter refers to local structure of objects. In differential geometry, a bridge between 
topology and geometry is given by the Gauss-Bonnet theorem. It states that for compact 
2D Riemannian manifolds M without boundary, the integral over the Gaussian curvature 
F(x) of the manifold is (i) integer and (ii) a topological invariant 

2(! -g) = y- [ d 2 xF(x). (8) 

2 vr J M 

Here, g is the genus of M, e.g., g — 0 for a 2D sphere and g = 1 for a 2D torus. The Gaussian 
curvature F(x) can be defined as follows. Attach to every point on M the tangential plane, 
a 2D vector space. Take some vector from the tangential plane at a given point on M and 
parallel transport it around an infinitesimal closed loop on M. The angle mismatch of the 
vector before and after the transport is proportional to the Gaussian curvature enclosed in 
the loop. 

In the physical systems that we want to describe, the manifold M is the BZ and the 
analogue of the tangent plane on M is a space spanned by the Bloch states of the occupied 
bands at a given momentum k E BZ. The Gaussian curvature of differential geometry is 
now generalized to a curvature form, called Berry curvature F. In our case, it is given by 
an n x n matrix of differential forms that is defined via the Berry connection A as 


F := -Fy(fc) d fcj A dkj (9a) 

F ij( k ) ■= d i A j{k) ^djA^k) + [^(fc),^*:)], i,j = 1,■••,<*, (9b) 

A:=A i (fc)dfc i , (9c) 

N 

Af > (k):=Y, U L(mU c , l ,(k), a,6=1,■■■,«, i = l,•■■,<!. (9d) 

a=1 
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(Two different conventions for the Berry connection are commonly used: Either it is purely 
real or purely imaginary. Here we choose the latter option.) The unitary transformation 
U(k) that diagonalizes the Hamiltonian was defined in Eq. (2b), both A i (k) and F-(k) 
are n x n matrices, we write d i = d/dk i and the sum over repeated spatial coordinate 
components i,j is implicit. 

Under a local U(n) gauge transformation in momentum space that acts on the states of 
the lower bands and is parametrized by the n x n matrix G(k) 

u aa( k ) — >U ab {k)G ba {k ), a — 1, - • ■ ,N, a = l,---,n, (10a) 

the Berry connection A changes as 

A —► G^AG + G'dG, (10b) 

while the Berry curvature F changes covariantly 

F —»■ G^FG, (10c) 


leaving its trace invariant. 

a. Chern numbers For the spatial dimension d = 2, the generalization of the Gauss- 
Bonnet theorem (8) in algebraic topology was found by Chern to be 

2C (1) :=— [ tr F 

2 nj 


BZ 


= 2— / d 2 fctrU 


2tt 


12 ' 


( 11 ) 


BZ 

This defines a gauge-invariant quantity, the first Chern number C (1) . Remarkably, C ( A can 
only take integer values. In order to obtain a topological invariant for any even dimension 
d = 2s of space, we can use the s-th power of the local Berry curvature form F (using the 
wedge product) to build a gauge invariant ri-form that can be integrated over the BZ to obtain 
scalar. Upon taking the trace, this scalar is invariant under the gauge transformation (10a) 
and defines the s-th Chern number 


2cW: 4Gl)7 trri ’ (i2) 

BZ 

where F s = F A ■ ■ ■ A F. As with the case s = 1 that we have exemplified above, is 
integer for any s = 1, 2, ■ • •. 
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From inspection of Tab. I we see that symmetry classes without chiral symmetry may 
have integer topological invariants Z only when the dimension d of space is even. In fact, all 
the integer invariants of these classes are given by the Chern number of the respective 
dimension. 

b. Winding numbers Let us now consider systems with chiral symmetry C. To con¬ 
struct their topological invariants as a natural extension of the above, we consider a different 
representation of the Chern numbers C (s b In terms of the flatland projector Hamiltonian 
Q(k) that was defined in Eq. (4a), we can write 

C ( *>oce h ... td Jd d ktT[Q(k)d tt Q(k)---d id Q(k)\, d = 2s. (13) 

BZ 

The form of Eq. (13) allows to interpret as the winding number of the unitary trans¬ 
formation Q(k) over the compact BZ. One verifies that O s ' = 0 for symmetry classes with 
chiral symmetry by inserting CC' at some point in the expression and anticommuting C 
with all Q, using the cyclicity of the trace. After 2s + 1 anticommutations, we are back 
to the original expression up to an overall minus sign and found = — C^. Hence, all 
systems with chiral symmetry have vanishing Chern numbers. 

In odd dimensions of space, we can define an alternative topological invariant for systems 
with chiral symmetry by modifying Eq. (13) and using the chiral operator C 

w ' s) := 2 ^ 1)1 d “ ktI [CQWS il Q(k)---d id a(k)} 

BZ 

= (2s + 1)! (^) £il - id / ddk tr ■ • • d id q{k)] , d = 2s + 1. 

BZ 

(14) 

Upon anticommuting the chiral operator C once with all matrices Q and using the cyclicity 
of the trace, one finds that the expression for W^ vanishes for even dimensions. The second 
line of Eq. (14) allows to interpret W (,s ' as the winding number of the unitary off-diagonal 
part q(k) of the chiral Hamiltonian that was defined in Eq. (7a). With Eq. (14) we have 
given topological invariants for all entries Z in odd dimensions d in Tab. I. 

In summary, we have now given explicit formulas for the topological invariants for all 
entries Z in Tab. I for systems with translational invariance. It is important to remember 
that the classification of Tab. I is restricted to systems without interactions. If interactions 
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are allowed, that neither spontaneously nor explicitly break the defining symmetry of a 
symmetry class, one of two things can happen: i) Two phases which are distinguished 
by a noninteracting invariant like W (0 ^ might, sometimes but not always, be connected 
adiabatically (i.e., without a closing of the spectral gap) by turning on strong interactions, 
ii) Interactions can enrich the classification of Tab. I by inducing new phases with topological 
response functions that are distinct from those of the noninteracting phases. We will given 
an example for the scenario i) in Sec IE. 

Besides, interactions can strongly modify the topological boundary modes of the nonin¬ 
teracting systems to the extend that they can be gapped without breaking the protective 
symmetries, but at the expense of introducing topological order on the boundary. 


C. The Su-Schrieffer-Heeger model 


The first example of a topological band insulator that we consider here is also the simplest: 
The Su-Schrieffer-Heeger model ’ describes a ID chain of atoms with one (spinless) electronic 
orbital each at half filling. The model was originally proposed to describe the electronic 
structure of poly acetylene. This ID organic molecule features a Peierls instability by which 
the hopping integral between consecutive sites is alternating between strong and weak. 
This enlarges the unit cell to contain two sites A and B. The second-quantized mean-field 
Hamiltonian reads 


N 

H = t [(! - s )cA,i c B,i + (1 + <S)4,i c x,i+i + h e. 


(15) 


2—1 

Here, c Ai and c^ Bi create an electron in the i-th unit cell on sublattice A and B , respec¬ 
tively. If we identify % — N + 1 = 1, periodic boundary conditions are implemented. The 
corresponding Bloch Hamiltonian 


H — t ^ Ca t kha/3,kCl3,k 

feeBZ a=A,B 


hk = 


0 (1 -6) + (1 + 6)e~ ik 

(1 - 5) + (1 + 5)e lk 0 

= a x [(1 — 6) + (1 + 5) cos k\ + a y (l + 5) sin k, 


(16a) 

(16b) 

(16c) 


where a x and a y are the first two Pauli matrices acting on the sublattice index, t is the 
nearest-neighbor hopping integral, and 5 is a dimensionless parametrization of the strong- 
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(a) 


(b) 



FIG. 1. Energy spectra for the Su-Schrieffer-Heeger model with open boundary conditions (a) 
in the trivial phase and (b) in the nontrivial topological phase with a zero energy mode on each 
boundary point. 

weak dimerization of bonds. 

We observe that Hamiltonian (16) has time-reversal symmetry T = /C, chiral symmetry 
C = u z and thus also particle-hole symmetry V = cr z IC. This places it in class BDI of Tab. I 
with a Z topological characterization. Observe that breaking the time-reversal symmetry 
would not alter the topological properties, as long as the chiral symmetry is intact. The 
model would then belong to class AIII, which also features a Z classification. Hence, it is the 
chiral symmetry that is crucial to protecting the topological properties of Hamiltonian (16). 
Notice that generic longer-range hopping (between sites of the same sublattice) breaks the 
chiral symmetry. 

What are the different topological sectors that can be accessed by tuning the parameter 
5 in the Su-Schrieffer-Heeger model? We observe that the dispersion 

e 2 k — 2 [(1 + 5 2 ) + (1 - 5 2 ) cos k] (17) 

is gapless for 5 = 0, hinting that this is the boundary between two distinct phases 5 > 0 and 
5 < 0. As we are interested in understanding the topological properties of these phases, we 
can analyze them for any convenient value of the parameter 6 and then conclude that they 
are the same in the entire phase by adiabaticity. We consider the Hamiltonian (15) with 
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open boundary conditions and choose the representative parameters 

• 5 = +1 : The operators c\ A and cjy B do not appear in the Hamiltonian for the open 
chain. Hence, there exists a state at either end of the open chain that can be occupied 
or unoccupied at no cost of energy. Thus, either end of the chain supports a localized 
topological end state (see Fig. 1). Away from 5 = +1, as long as 5 > 0, the end 
states start to overlap and split apart in energy by an amount that is exponentially 
small in the length N of the chain. We can back np this observation by evaluating the 
topological invariant (14) for this phase. The off-diagonal projector is qk = e~ lk and 
its winding number evaluates to 

W (0) = — f dk e ifc (—i)e -ifc = 1. (18) 

2tt J 

• 6 = — 1 : In this case strong bonds form between the two sites in every unit cell and no 
topological end states appear. Correspondingly, as the off-diagonal projector qk = 1 
is independent of k, we conclude that the winding number vanishes W (0) = 0. 

One can visualize the winding number of a two-band Hamiltonian that has the form 
hk — dk ■ cr in the following way. If the Hamiltonian has chiral symmetry, we can choose 
this symmetry to be represented by C = u z without loss of generality. Then dk has to lie in 
the x-y-plane for every k and may not be zero if the phase is gapped. The winding number 
W !0) measures how often dk winds around the origin in the x-y-plane as k changes from 0 
to 2ir. 

Besides the topological end states, the Su-Schrieffer-Heeger model also features topolog¬ 
ical domain wall states between a region with 5 > 0 and S < 0. Such topological midgap 
modes have to appear pairwise in any periodic geometry. As the system is considered at 
half filling, each of these modes binds half an electron charge. This is an example of charge 
fractionalization at topological defects. It is important to remember that these defects are 
not dynamical, but are rigidly fixed external perturbations. Therefore, this form of fraction- 
alization is not related to intrinsic topological order. 

D. The one-dimensional p- wave superconductor 

In the Su-Schrieffer-Heeger model, particle-hole symmetry (and with it the chiral symme¬ 
try) is in some sense fine-tuned, as it is lost if generic longer-range hoppings are considered. 
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In superconductors, particle-hole symmetry arises more naturally as a symmetry that is 
inherent in the redundant description of mean-held Bogoliubov-deGennes Hamiltonians. 

Here, we want to consider the simplest model for a topological superconductor that has 
been studied by Kitaev in Ref. 6. The setup is again a ID chain with one orbital for spinless 
fermion on each site. Superconductivity is encoded in pairing terms c\c\ +1 that do not 
conserve particle number. The Hamiltonian is given by 

N 


H = 


Y ( C i c i+1 + 4+iCi) - jCi + M+d + A‘CiC i+1 


(19) 


i= 1 

Here, /i is the chemical potential and A is the superconducting order parameter, which we 
will decompose into its amplitude |A| and complex phase 9, i.e., A = |A|e 10 . 

The fermionic operators c\ obey the algebra 


{ C i ) C j 1 


( 20 ) 


with all other anticommutators vanishing. We can chose to trade the operators c\ and on 
every site i for two other operators a.; and fe; that are defined by 


a ?; = e ie/2 Ci + e l6/2 c\, b t = y (e l6/2 Ci - e 10/2 cfj . 
These so-called Majorana operators obey the algebra 


( 21 ) 


{ai,aj} = {bi,b :j } = 2%, 1.) = (} VtJ. 


( 22 ) 


In particular, they square to 1 



(23) 


and are self-conjugate 

a\ = di, b\ — bi. (24) 

In fact, we can always break up a complex fermion operator on a lattice site into its real 
and imaginary Majorana components though it may not always be a useful representation. 
As an aside, note that the Majorana anti-commutation relation in Eq. (22) is the same as 
that of the generators of a Clifford algebra where the generators all square to +1. Thus, 
mathematically one can think of the operators a* (or bi) as matrices forming by themselves 
the representation of Clifford algebra generators. 
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FIG. 2. Schematic illustration of the lattice p-wave superconductor Hamiltonian in the (a) trivial 
limit (b) non-trivial limit. The white (empty) and red (filled) circles represent the Majorana 
fermions making up each physical site (oval). The fermion operator on each physical site ( Cj ) is 
split up into two Majorana operators (a ? - and bj). In the non-trivial phase the unpaired Majorana 
fermion states at the end of the chain are labelled with a\ and bjy. These are the states which are 
continuously connected to the zero-modes in the non-trivial topological superconductor phase. 


When rewritten in the Majorana operators, Hamiltonian (19) takes (up to a constant) 
the form 


. N 

H — - [— iia>i bi + (t + | A|)6* a i+ i + (~t + | A|)a ? : &*+ 1 ]. (25) 

i=l 

After imposing periodic boundary conditions, it is again convenient to study the system in 
momentum space. When defining the Fourier transform of the Majorana operators a; = 
'Yhi elkla k we note that the the self-conjugate property (24) that is local in position space 
translates into a' k = a- k in momentum space (and likewise for the b k )- The momentum 
space representation of the Hamiltonian is 


H 


y y (a k b k )h k 

/cGBZ a=A,B 



h k = 


0 


^ — it cos k + |A| sin k 




+ it cos k + |A| sin k 


= a^lAI sin k + a v ( — — t cos k 


(26a) 

(26b) 

(26c) 


While this Bloch Hamiltonian is formally very similar to that of the Su-Schrieffer-Heeger 
model (16), we have to keep in mind that it acts on entirely different single-particle degrees of 
freedom, namely in the space of Majorana operators instead of complex fermionic operators. 
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As with the case of the Su-Schrieffer-Heeger model, the Hamiltonian (26) has a time-reversal 
symmetry T = <J Z K. and a particle-hole symmetry P — JC which combine to the chiral 
symmetry C = a z . Hence, it belongs to symmetry class BDI as well. For the topological 
properties that we explore below, only the particle-hole symmetry is crucial. If time-reversal 
symmetry is broken, the model changes to symmetry class D, which still supports a Z 2 
topological grading. 

To determine its topological phases, we notice that Hamiltonian (26) is gapped except 
for \t\ = |/i/21. We specialize again on convenient parameter values on either side of this 
potential topological phase transition 

• fi — 0, |A| — t : The Bloch matrix takes exactly the same form as that of the 
Su-Schrieffer-Heeger model (16) for the parameter choice 6 = +1. We conclude that 
the Hamiltonian (26) is in a topological phase. The Hamiltonian reduces to 



(27) 


A pictorial representation of this Hamiltonian is shown in Fig. 2 b). With open 
boundary conditions it is clear that the Majorana operators a± and fejy are not coupled 
to the rest of the chain and are ‘unpaired’. In this limit the existence of two Majorana 
zero modes localized on the ends of the chain is manifest. 

• A — t — 0, n < 0 : This is the topologically trivial phase. The Hamiltonian is 
independent of k and we conclude that the winding number vanishes W (0) = 0. I 11 this 
case the Hamiltonian reduces to 



(28) 


In its ground state the Majorana operators on each physical site are coupled but the 
Majorana operators between each physical site are decoupled. In terms of the physical 


complex fermions, it is the ground state with either all sites occupied or all sites empty. 


A representation of this Hamiltonian is shown in Fig. 2 a). The Hamiltonian in the 
physical-site basis is in the atomic limit, which is another way to see that the ground 
state is trivial. If the chain has open boundary conditions there will be no low-energy 
states on the end of the chain if the boundaries are cut between physical sites. That 
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is, we are not allowed to pick boundary conditions where a physical complex fermionic 
site is cut in half. 

These two limits give the simplest representations of the trivial and non-trivial phases. 
By tuning away from these limits the Hamiltonian will have some mixture of couplings 
between Majorana operators on the same physical site, and operators between physical 
sites. However, since the two Majorana modes are localized at different ends of a gapped 
chain, the coupling between them will be exponentially small in the length of the wire and 
they will remain at zero energy. In fact, in the non-trivial phase the zero modes will not be 
destroyed until the bulk gap closes at a critical point. 

It is important to note that these zero modes count to a different many-body ground state 
degeneracy than the end modes of the Su-Schrieffer-Heeger model. The difference is rooted 
in the fact that one cannot build a fermionic Fock space out of an odd number of Majorana 
modes, because they are linear combinations of particles and holes. Rather, we can define 
a single fermionic operator out of both Majorana end modes cq and b^ as c' := a i + i bjy. 
The Hilbert space we can build out of a\ and bjy is hence inherently nonlocal. This nonlocal 
state can be either occupied or empty giving rise to a two-fold degenerate ground state of 
the chain with two open ends. (In contrast, the topological Su-Schrieffer-Heeger chain has 
a four-fold degenerate ground state with two open ends, because it has one fermionic mode 
on each end.) The Majorana chain thus displays a different form of fractionalization than 
the Su-Schrieffer-Heeger chain. For the latter, we observed that the topological end modes 
carry fractional charge. In the Majorana chain, the end modes are a fractionalization of a 
fermionic mode into a superposition of particle and hole (and have no well defined charge 
anymore), but the states |0) (with c|0) = 0) and <dj0) do have distinct fermion parity. The 
nonlocal fermionic mode formed by two Majorana end modes is envisioned to work as a 
qubit (a quantum-mechanical two-level system) that stores quantum information (its state) 
in a way that is protected against local noise and decoherence. 

E. Reduction of the 10-fold way classification by interactions: Z —> Zg in class BDI 

When time-reversal symmetry T = /C is present, the model considered in Sec. ID belongs 
to class BDI of the classification of noninteracting fermionic Hamiltonians in Tab. I with a 
Z topological characterization. We want to explore how interactions alter this classification, 
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following a calculation by Fidkowski and Kitaev from Ref. 8. To this end, we consider a 
collection of n identical ID topological Majorana chains in class BDI and only consider their 
Majorana end modes on one end, which we denote by 7 1; • • • , y n . We will take the point of 
view that if we can gap the edge, we can continue the bulk to a trivial state (insulator). This 
is not entirely a correct point of view in general (see 2D topologically ordered states such 
as the toric code discussed in the next Section), but works for our purposes. Given some 
integer n, we ask whether we can couple the Majorana modes locally on one end such that 
no gapless degrees of freedom are left on that end and the ground state with open boundary 
conditions becomes singly degenerate. To remain in class BDI, we only allow couplings that 
respect time-reversal symmetry. Let us first derive the action of T on the Majorana modes. 
The complex fermion operators are left invariant under time-reversal Tc'T 1 = c. Hence, 

T(a + ib)T= TaT- 1 - \TbT~ 1 = a + ib =7 T aT~ l = a, TbT~ x = -b. (29) 

Thus, when acting on the modes localized on the left end of the wire (which transform like 
the a’s), time-reversal symmetry leaves the Majorana operators invariant. 

The most naive coupling term that would gap out two Majoranas is iy^. This is because 
two Majoranas can form a local Hilbert space (unlike just one Majorana), and this local 
Hilbert space can be split unless some other symmetry prevents it from being split. However, 
time-reversal symmetry forbids these hybridization terms, for it sends vy\^ 2 —> —i7i7 2 - In 
spinful systems, another symmetry which can do this is MT, where M is a mirror operator 
(which in spinful systems squares to —1 M 2 = — 1) and T is the usual time-reversal operator 
T 2 = —1, such that {MT) 2 = M 2 T 2 = 1 and hence MT acts like spinless time reversal.' 
Realizing that such a term is not allowed is the end of the story for noninteracting systems, 
yielding the classification Z. Lets find out what interactions do to this system. The steps 
that we will now outline are summarized in Fig. 3. 

We saw that two Majorana end states cannot be gapped: the only possible interacting 
or noninteracting Hamiltonian is ia\a 2 . Three Majoranas clearly cannot be gapped either, 
as it is an odd number. Let us thus add two more Majorana end states into the mix. Any 
one-body term still is disallowed but the term 

Hint — aia 2 a 3 a 4 (30) 

can be present. We can now form two complex fermions, C\ = (cq + ia 2 )/\/2, c 2 = (a 3 + 
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2 wires 


tfint 


4 wires 



8 wires 


FIG. 3. Schematic illustration of the many body energy levels for 2, 4, and 8 wires with Majorana 
end states as well as the (partial) lifting of their degeneracy by the Hamiltonians in Eqs. (31) 
and (32). 


\a±)/\/2. In terms of these two fermions, the Hamiltonian reads 



where ri\ = c\c\ and n 2 = 4 C 2 are the occupation numbers. The Hamiltonian is diagonal 
in the eigenbasis \nin 2 ) of the occupation number operators, and the states 111), 100) are 
degenerate at energy —1/4, while the states |01), 110) are degenerate at energy +1/4. The 
original noninteracting system of four Majorana fermions had a degeneracy of 2 2 = 4. The 
interaction, however, has lifted this degeneracy, but not all the way to a single nondegen¬ 
erate ground state. Irrespective of the sign of the interaction, it leaves the states doubly 
degenerate on one edge, and hence cannot be adiabatically continued to the trivial state of 
single degeneracy. However, if we add four more Majoranas wires so that we have n — 8 
Majoranas, we can build an interaction which creates a singly degenerate ground state. We 
can understand this as follows: Add two interactions 


//j n l — —Qf(ciifl2®3®4 + 


(32a) 


These create two doublets, one in c\, c 2 defined above, and one in C 3 = ( 0,5 + ia e )/\/2, c 4 = 
( 0,7 + ias)/\/2. We couple these doublets via the interaction 


H int = J2 PiA 


Cl 


(4 cj)<7i 


C 3 


(32b) 


i=x,y,z \ / \ £4 J 

Representing each of the doublets as a spin-1/2 S, this interaction is nothing but an S ■ S 
term. If we take 0 < /3 <C a, then we can approximate the interaction (3 by its action within 
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the two ground state doublets. As such, this interaction creates a singlet and a triplet (in 
that doublet) and for the right sign of f3, we can put the singlet below the triplet, thereby 
creating a unique ground state 

-)=(|0110)-11001)), (33) 

in terms of the occupation number states |nin 2 n 3 n 4 ). This unique ground state can be 
adiabatically continued to the atomic limit. In this way the noninteracting Z classification 
of class BDI breaks down to Z§ if interactions are allowed. 
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II. EXAMPLES OF TOPOLOGICAL ORDER 


So far, we have been concerned with symmetry protected topological states and consid¬ 
ered examples that were motivated by the topological classification of free fermion Hamilto¬ 
nians. The topological properties of these systems are manifest by the presence of protected 
boundary modes. 

In this Section, we want to familiarize ourselves with the concept of intrinsic topological 
order by ways of two examples. We will study the connections between different characteri¬ 
zations of topological order, such as fractionalized excitations in the bulk and the topological 
ground state degeneracy. Our examples will be in 2D space, as topologically ordered states 
do not exist in ID and are best understood in 2D. Our first example, the toric code, has 
Abelian anyon excitations, while the second example, the chiral p -wave superconductor, 
features non-Abelian anyons. 


A. The toric code 

The first example of a topologically ordered state is an exactly soluble model with van¬ 
ishing correlation length. The significance of having zero correlation length is the following. 
The correlation functions of local operators decay exponentially in gapped quantum ground 
states in ID and 2D with a characteristic length scale given by the correlation length £. 9 
In contrast, topological properties are encoded in quantized expectation values of nonlocal 
operators (for example the Hall conductivity) or the degeneracy of energy levels (such as 
the end states of the Su-Schrieffer-Heeger model). In finite systems, such quantizations and 
degeneracies are generically only exact up to corrections that are of order e ~ L ^, where L 
is the linear system size. Models with zero correlation length are free from such exponen¬ 
tial finite-size corrections and thus expose the topological features already for the smallest 
possible system sizes. The down-side is that their Hamiltonians are rather contrived. 

We define the toric code model 10 on a square lattice with a spin-1/2 degree of freedom on 
every bond j (see Fig. 4). The four spins that sit on the bonds emanating from a given site 
of the lattice are referred to as a star s. The four spins that sit on the bonds surrounding a 
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FIG. 4. The toric code model is defined on a square lattice with spin-1/2 degrees of freedom on 
every bond (black squares). The operator A s acts with a x on all four spins one the bonds that 
are connected to a lattice site (a star s). The operator B p acts with a z on all four spins around a 
plaquette p. 


square of the lattice are called a plaquette p. We define two sets of operators 



(34) 


that act on the spins of a given star s and plaquette p, respectively. Here, cr x,z are the 
respective Pauli matrices acting on the spin on bond j. 

These operators have two crucial properties which are often used to construct exactly 
soluble models for topological states of matter 

1. All of the A s and B p commute with each other. This is trivial for all cases except 
for the commutator of A s with B p if s and p have spins in common. However, any 
star shares with any plaquette an even number of spins (edges), so that commuting 
A s with B p involves commuting an even number of a z with a x , each of which comes 
with a minus sign. 

2. The operators 


1 ~B P 


(35) 


2 


2 


are projectors. The former projects out plaquette states with an even number of spins 
polarized in the positive ^-direction. The latter projects out stars with an even number 
of spins in the positive x-direction. 
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1. Ground states 


The Hamiltonian is defined as a sum over these commuting projectors 

H=-J e J2A s - ( 36 ) 

s p 

where the sums run over all stars s and plaquettes p of the lattice. Let us assume that both 
J e and J m are positive constants. Then, the ground state is given by a state in which all 
stars s and plaquettes p are in an eigenstate with eigenvalue +1 of A s and B p , respectively. 
(The fact that all A s and B p commute allows for such a state to exist, as we can diagonalize 
each of them separately.) Let us think about the ground state in the eigenbasis of the a x 
operators and represent by bold lines those bonds with spin up and and draw no lines along 
bonds with spin down. Then, A s imposes on all spin configurations with nonzero amplitude 
in the ground state the constraint that an even number of bold lines meets at the star s. 
In other words, we can think of the bold lines as connected across the lattice and they may 
only form closed loops. Bold lines that end at some star (“open strings”) are not allowed 
in the ground state configurations; they are excited states. Having found out which spin 
configurations are allowed in the ground state, we need to determine their amplitudes. This 
can be inferred from the action of the B p operators on these closed loop configurations. The 
B p flips all bonds around the plaquette p. Since B 2 = 1, given a spin configuration |c) in 
the cr^-basis, we can write an eigenstate of B p with eigenvalue 1 as 

-j=(\c) + B p \c)), (37) 

for some fixed p. This reasoning can be extended to all plaquettes so that we can write for 
the ground state 



where |c) is a closed loop configuration [see Fig. 5 a)]. Is |GS) independent of the choice of 
|c>? In other words, in the ground state unique? We will see that the answer depends on 
the topological properties of the manifold on which the lattice is defined and thus reveals 
the topological order imprinted in |GS). 

To answer these questions, let us consider the system on two topologically distinct mani¬ 
folds, the torus and the sphere. To obtain a torus, we consider a square lattice with L x x L y 



+ ••• 


a) 

|GS) = 








□ 


b) 




FIG. 5. Visualization of the toric code ground states on the torus, a) The toric code ground state is 
the equal amplitude superposition of all closed loop configurations. b)-e) Four base configurations 
|c) entering Eq. (38) that yield topologically distinct ground states on the torus. 


sites and impose periodic boundary conditions. This lattice hosts 2L x L y spins (2 per unit 
cell for they are centered along the bonds). Thus, the Hilbert space of the model has di¬ 
mension 2 2LxLy . There are L x L y operators A s and just as many B p . Hence, together they 
impose 2 L x L y constraints on the ground state in this Hilbert space. However, not all of 
these constraints are independent. The relations 


i=n*. i= n^ 


(39) 


make two of the constraints redundant, yielding (2 L x L y — 2) independent constraints. The 
ground state degeneracy (GSD) is obtained as the quotient of the Hilbert space dimension 
and the subspace modded out by the constraints 


~)‘2L x I j u 


GSD = 


C 2^LxLy 2 


= 4. 


(40) 


The four ground states on the torus are distinguished by having an even or an odd number 
of loops wrapping the torus in the x and y direction, respectively. Four configurations |c) 
that can be used to build the four degenerate ground states are shown in Fig. 5 b)-e). 
This constitutes a set of “topologically degenerate” ground states and is a hallmark of the 
topological order in the model. 

Let us contrast this with the ground state degeneracy on the sphere. Since we use a 
zero correlation length model, we might as well use the smallest convenient lattice with the 
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topology of a sphere. We consider the model (36) defined on the edges of a cube. The 
same counting as above yields that there are 12 degrees of freedom (the spins on the 12 
edges), 8 constraints from the A s operators defined on the corners and 6 constraints from 
the B p operators defined on the faces. Subtracting the 2 redundant constraints (39) yields 
12 — (8 + 6 — 2) = 0 remaining degrees of freedom. Hence, the model has a unique ground 
state on the sphere. 

On a general manifold, we have 


GSD 


2number of noncontractible loops 


(41) 


An important property of the topologically degenerate ground states is that any local oper¬ 
ator has vanishing off-diagonal matrix elements between them in the thermodynamic limit. 
Similarly, no local operator can be used to distinguish between the ground states. We can, 
however, define nonlocal operators that transform one topologically degenerate ground state 
into another and that distinguish the ground states by topological quantum numbers. (No¬ 
tice that such operators may not appear in any physical Hamiltonian due to their nonlocality 
and hence the degeneracy of the ground states is protected.) On the torus, we define two 
pairs of so-called Wilson loop operators as 

'UV := n V 1 . W “/y n <*• ( 42 > 

i&k/v i&x/y 

Here, l e x , are the sets of spins on bonds parallel to a straight line wrapping the torus 
once along the x- and y-direction, respectively. The l™/ y are the sets of spins on bonds 
perpendicular to a straight line that connects the centers of plaquettes and wraps the torus 
once along the x and y-direction, respectively. We note that the ITT and W? /y commute 
with all A s and B p 



and thus also with the Hamiltonian. Furthermore, they obey 

W e x w™ = -w™w°. ( 44 ) 

This algebra must be realized in any eigenspace of the Hamiltonian. However, due to 
Eq. (44), it cannot be realized in a one-dimensional subspace. We conclude that all 
eigenspaces of the Hamiltonian, including the ground state, must be degenerate. In the 
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FIG. 6. Visualization of operations to compute the braiding statistics of toric code anyons. a) 
Two e excitations above the ground state, b) Two m excitations above the ground state, c) Loop 
created by braiding ei around e, 2 . c) Loop created by braiding e around m. A phase of —1 results 
for this process because there is a single bond on which both a a x operator (dotted line) and a a z 
operator (bold line) act. 

a x basis that we used above, W™ y measures whether the number of loops wrapping the 
torus is even or odd in the x and y direction, respectively, giving 4 degenerate ground states. 
In contrast, W£/ y changes the number of loops wrapping the torus in the x and y direction 
between even and odd. 

2. Topological excitations 

To find the topological excitations of the system above the ground state, we ask which 
are the lowest energy excitations that we can build. Excitations are a violation of the rule 
that all stars s are eigenstates of A s and all plaquettes p are eigenstates of B p . Let us first 
focus on star excitations which we will call e. They appear as the end point of open strings, 
i.e., if the closed loop condition is violated. Since any string has two end points, the lowest 
excitation of this type is a pair of e. They can be created by acting on the ground state 
with the operator 

w?:=n«J. < 45 ) 

je/ e 

where l e is a string of bonds connecting the two excitations e\ and e 2 [see Fig. 6 a)]. The 
state 

|er,e 2 ) := W£|GS) (46) 
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has energy 4 J e above the ground state energy. Similarly, we can define an operator 

:= n y. < 47 ) 

that creates a pair of plaquette defects mi and m 2 connected by the string l m of perpendicular 
bonds [see Fig. 6 b)]. (Notice that the operator W t ™ does not flip spins when the ground state 
is written in the a x basis. Rather, it gives weight +1/ — 1 to the different loop configurations 
in the ground state, depending on whether an even or an odd number of loops crosses l m .) 
The state 

| mi , m2 ) : = W”|GS> (48) 

has energy 4J m above the ground state energy. Notice that the excited states |ei,e 2 ) and 
|m!,m 2 ) only depend on the positions of the excitations and not on the particular choice of 
string that connects them. Furthermore, the energy of the excited state is independent of 
the separation between the excitations. The excitations are thus “deconfined”, i.e., free to 
move independent of each other. 

It is also possible to create a combined defect when a plaquette hosts a m excitation and 
one of its corners hosts a e excitation. We call this combined defect / and formalize the 
relation between these defects in a so-called fusion rule 

e x m = /. (49a) 

When two e-type excitations are moved to the same star, the loop l e that connects them 
becomes a closed loop and the state returns to the ground state. For this, we write the 
fusion rule 

e x e = 1, (49b) 

where 1 stands for the ground state or vacuum. Similarly, moving two m-type excitations 
to the same plaquette creates a closed loop l m , which can be absorbed in the ground state, 
i.e., 

m x m — 1. (49c) 

Superimposing the above processes yields the remaining fusion rules 

m x / = e, e x / = m, f x f — 1. (49d) 

It is now imperative to ask what type of quantum statistics these emergent excitations 
obey. We recall that quantum statistics are defined as the phase by which a state changes if 
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two identical particles are exchanged. Rendering the exchange operation as an adiabatically 
slow evolution of the state, in three and higher dimensions only two types of statistics are 
allowed between point particles: that of bosons with phase +1 and that of fermions with 
phase —1. In 2D, richer possibilities exist and the exchange phase 6 can be any complex 
number on the unit circle, opening the way for any ons. While the exchange is only defined for 
quantum particles of the same type, the double exchange (braiding) is well defined between 
any two deconhned anyons. We can compute the braiding phases of the anyons e, m, and / 
that appear in the toric code one by one. Let us start with the phase resulting from braiding 
e\ with e 2 . The initial state is IF;®|GS) depicted in Fig. 6 a). Moving e\ around e 2 leaves a 
loop of flipped <j x bonds around e 2 [see Fig. 6 c)]. This loop is created by applying B p to all 
plaquettes enclosed by the loop l* along which e\ moves. We can thus write the final state 
as 



W, e .|GS) 


IF® 



| GS) 


= IF®|GS). 


(50) 


Flipping the spins in a closed loop does not alter the ground state as it is the equal amplitude 
of all loop configurations. We conclude that the braiding of two e particles gives no phase. 
Similar considerations can be used to conclude that the braiding of two m particles is trivial 
as well. In fact, not only the braiding, but also the exchange of two e particles and two m 
particles is trivial. (We have not shown that here.) 

More interesting is the braiding of m with e. Let the initial state be IF™Wq®[CIS) and 
move the e particle located on one end of the string /® n around the magnetic particle m on 
one end of the string l m . Again this is equivalent to applying B p to all plaquettes enclosed 
by the path Z® of the e particle, so that the final state is given by 



IF^IFfe|GS> 


m 



IFfe|GS) 


- IF“IF®|GS>. 


(51) 


The product over B p operators anticommutes with the path operator IF/S, because there 
is a single bond on which a single a x and a single <j z act at the crossing of Z m and Z® [see 
Fig. 6 d)]. As a result, the initial and final state differ by a —1, which is the braiding phase 
of e with m. Particles with this braiding phase are called (mutual) semions. 
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Notice that we have moved the particles on contractible loops only. If we create a pair of e 
or m particles, move one of them along a noncontractible loop on the torus, and annihilate the 
pair, we have effectively applied the operators and W? /y to the ground state (although 
in the process we have created finite energy states). The operation of moving anyons on 
noncontractible loops thus allows to operate on the manifold of topologically degenerate 
groundstates. This exposes the intimate connection between the presence of fractionalized 
excitations and topological groundstate degeneracy in topologically ordered systems. 

From the braiding relations of e and m we can also conclude the braiding and exchange 
relations of the composite particle /. This is most easily done in a pictorial way by repre¬ 
senting the particle worldlines as moving upwards. For example, we represent the braiding 
relations of e and m as 



e e e e 

The exchange of two /, each of which 



e m e m 
is composed of one e and one m is then 



jl n 

me 


Ji n 

me 


(52) 


(53) 


Notice that we have used Eq. (52) to manipulate the crossing in the dotted rectangles. 
Exchange of two / thus gives a phase —1 and we conclude that / is a fermion. 

In summary, we have used the toric code model to illustrate topological ground state 
degeneracy and emergent anyonic quasiparticles as hallmarks of topological order. We note 
that the toric code model does not support topologically protected edge states. 


B. The two-dimensional p-wave superconductor 

The second example of a 2D system with anyonic excitations that we want to discuss 
here is the chiral p-wave superconductor. Unlike the toric code, due to its chiral nature, it is 
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a model with nonzero correlation length. The vortices of the chiral p -wave superconductor 
exhibit anyon excitations which have exotic non-Abelian statistics. 12-14 (The anyons in the 
toric code are Abelian, we will see below what that distinction refers to.) For the system to 
be topologically ordered, these vortices should appear as emergent, dynamical excitations. 
This requires to treat the electromagnetic gauge field quantum-mechanically. (In fact, since 
the fermion number conservation is spontaneously broken down to the conservation of the 
fermion parity in the superconductor, the relevant gauge theory involves only a Z 2 instead 
of a U(l) gauge field.) However, the topological properties that we want to discuss here 
can also be seen if we model the gauge field and vortices as static defects, rather than 
within a fluctuating Z 2 gauge theory. This allows us to study a models very similar to 
the “noninteracting” topological superconductor in ID and still expose the non-Abelian 
statistics. 

For pedagogy we will use both lattice and continuum models of the chiral superconductor. 
We begin with the lattice Hamiltonian defined on a square lattice 


H = 




t ( Dn+l.nDn,! 


n ' r* 

1 m,n+lm,n 


+ h.c. j — (n — 4 t)c. 


r 

m,nm,n 


(54) 


+ ( Ac t + 1,„4.,» + iAc L,„+l C ln,n + hx.)| ■ 

The fermion operators c m:7l annihilate fermions on the lattice site (m, n) and we are consid¬ 
ering spinless (or equivalently spin-polarized) fermions. We set the lattice constant a = 1 
for simplicity. The pairing amplitude is anisotropic and has an additional phase of i in the 
^/-direction compared to the pairing in the x-direction. Because the pairing is not on-site, 
just as in the lattice version of the p -wave wire, the pairing terms will have momentum de¬ 
pendence. We can write this Hamiltonian in the Bogoliubov-deGennes form and, assuming 
that A is translationally invariant, can Fourier transform the lattice model to get 


H BdG — ^ 

p 

where e(p) = — 2f(cosp x + cos p y ) — (ji — 4 1) and = ^c p . For convenience we have 

shifted the chemical potential by the constant 4 1. As a quick aside we note that the model 
takes a simple familiar form in the continuum limit (p —>• 0): 


e(p) 2iA(sinp a; + isinpy) 

-2iA*(sinp3, — i sinp y ) ~ e (p) 


4b 


(55) 


H 


(cont) _ 

BdG — 



~ I 1 2iA (P* + i Py) \ x]) 

-2iA* (p x -ip y ) -^ + // j 


(56) 
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where m = 1/2 1 and p 2 = p 2 x + p 2 . We see that the continuum limit has the characteristic 
Px + i P y chiral form for the pairing potential. The quasiparticle spectrum of is 

E± = ±a/ (p 2 /2m — p) 2 + 4| A| 2 p 2 , which, with a nonvanishing pairing amplitude, is gapped 
across the entire BZ as long as p 7 ^ 0. This is unlike some other types of p -wave pairing 
terms [e.g., A (p) = A p x \ which can have gapless nodal points or lines in the BZ for p > 0. 
In fact, nodal superconductors, having gapless quasiparticle spectra, are not topological 
superconductors by definition (i.e., a bulk excitation gap does not exist). 

We recognize the form of as a massive 2D Dirac Hamiltonian, and indeed Eq. (54) 

is just a lattice Dirac Hamiltonian which is what we will consider first. In the first quantized 
notation, the single particle Hamiltonian for a superconductor is equivalent to that of an 
insulator with an additional particle-hole symmetry. It is thus placed in class D of Tab. I 
and admits a Z topological classification in 2D. Thus, we can classify the eigenstates of 
Hamiltonian (54) by a Chern number - but due to the breaking of U(l) symmetry, the 
Chern number does not have the interpretation of Hall conductance. However, it is still a 
topological invariant. 

We expect that Hb<ig will exhibit several phases as a function of A and p for a fixed t > 0. 
For simplicity let us set t = 1/2 and make a gauge transformation c p —> e l6 ^ 2 c p , cj, —> e~ ie,/2 c t 
where A = \A\e l6 . The Bloch Hamiltonian for the lattice superconductor is then 

■HbcIgCp) = (2 - p - cos p x - cos p y ) < 7 z - 2|A| sin p x a y - 2| A| sin p y a x , (57) 

where the 07 , i = x, 7 /,z, are the Pauli matrices in the particle/hole basis. Assuming 
|A| 7 ^ 0, this Hamiltonian has several fully-gapped superconducting phases separated by 
gapless critical points. The quasi-particle spectrum for the lattice model is 

E± — ± \J(2 — p — cos p x — cos p y ) 2 + 4| A | 2 sin 2 ^ + 4|A | 2 sin 2 p y (58) 

and is gapped (under the assumption that |A| 7 ^ 0) unless the prefactors of all three Pauli 
matrices vanish simultaneously. As a function of (p x ,Py,p) we End three critical points. 
The first critical point occurs at (p x ,P y ,p) = (0,0,0). The second critical point has two 
gap-closings in the BZ for the same value of p : (n, 0, 2) and (0, 7 r, 2). The third critical 
point is again a singly degenerate point at ( 7 r, 7r, 4). We will show that the phases for p < 0 
and p > 4 are trivial superconductors while the phases 0 < p < 2 and 2 < p < 4 are 
topological superconductors with opposite chirality. In principle one can define a Chern 
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number topological invariant constructed from the eigenstates of the lower quasi-particle 
band to characterize the phases. We will show this calculation below, but first we make 
some physical arguments as to the nature of the phases, following the discussion in Ref. 11. 

We will first consider the phase transition at // = 0. The low-energy physics for this 
transition occurs around ( p x ,Py ) = (0,0) and so we can expand the lattice Hamiltonian 
around this point; this is nothing but Eq. (56). One way to test the character of the /i < 0 
and /i > 0 phases is to make an interface between them. If we can find a continuous 
interpolation between these two regimes which is always gapped then they are topologically 
equivalent phases of matter. If we cannot find such a continuously gapped interpolation 
then they are topologically distinct. A simple geometry to study is a domain wall where 
/i = n(x) such that /j,(x ) = —po for x < 0 and p(x) = +po for x > 0 for a positive constant 
/io. This is an interface which is translationally invariant along the ^/-direction, and thus we 
can consider the momentum p y as a good quantum number to simplify the calculation. What 
we will now show is that there exist gapless, propagating fermions bound to the interface 
which prevent us from continuously connecting the p < 0 phase to the p > 0 phase. This is 
one indication that the two phases represent topologically distinct classes. 

The single-particle Hamiltonian in this geometry is 


^BdG (p y ) - ^ 


-H(x) 

-2i|A| (~i£ 


1 Py) 


2i|A| 


+ [ Py) 

H(x) 


(59) 


where we have ignored the quadratic terms in p, and p y is a constant parameter, not an op¬ 
erator. This is a quasi-ID Hamiltonian that can be solved for each value of p y independently. 
We propose an ansatz for the gapless interface states: 


I ^py(x,y)) = e iPyV exp ^-^j- p(x')dx'^J \(j> 0 ) (60) 

for a constant, normalized spinor |0o) - The secular equation for a zero-energy mode at p y = 0 

is 


'HBdG(Py)H’o( x ^y)) = 0 


—p(x) —fj>(x) 
n(x) p(x) 


I0o> = o. 


(61) 


The constant spinor which is a solution of this equation is |0o) = l/\/2 (1, — 1) T . This form 
of the constant spinor immediately simplifies the solution of the problem at finite p y . We see 
that the term proportional to p y in Eq. (59) is —2|A| p y a x . Since cr^l^o) = — |A>o)- he-, the 
solution |0 O ) is an eigenstate of a x , we conclude that \ijj py (x, y)) is an eigenstate of 'HBdG(/b/) 
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with energy E(p y ) = —2|A| p y . Thus, we have found a normalizable bound state solution at 
the interface of two regions with p < 0 and p > 0 respectively. This set of bound states, 
parameterized by the conserved quantum number p y is gapless and chiral, i.e., the group 
velocity of the quasiparticle dispersion is always negative and never changes sign (in this 
simplified model). The chirality is determined by the sign of the “spectral” Chern number 
mentioned above which we will calculate below. 

These gapless edge states have quite remarkable properties and are not the same chiral 
complex fermions that propagate on the edge of integer quantum Hall states, but chiral real 
(Majorana) fermions. Using Clifford algebra representation theory it can be shown that 
the so-called chiral Majorana (or Majorana-Weyl) fermions can only be found in spacetime 
dimensions (8k + 2), where k = 0,1, 2, • • ■ . Thus, we can only find chiral-Majorana states in 
(1 + 1) dimensions or in (9 + 1) dimensions (or higher!). In condensed matter, we are stuck 
with (1 + 1) dimensions where we have now seen that they appear as the boundary states 
of chiral topological superconductors. The simplest interpretation of such chiral Majorana 
fermions is as half of a conventional chiral fermion, i.e., its real or imaginary part. To show 
this, we will consider the edge state of a Chern number 1 quantum Hall system for a single 
edge 


= hv J2 pr fl^ ( 62 ) 

p 

where p is the momentum along the edge. The fermion operators satisfy jr/j,, r/ p , j = 5 pp >. 
Similar to the discussion on the ID superconducting wire we can decompose these operators 
into their real and imaginary Majorana parts 


Vp = 2 (7i ,p + 172 , P ), vl = 2 (7i,-p - 172 ,- P ), 


(63) 


where 7 0)P (a = 1, 2) are Majorana fermion operators satisfying 7 ^ = 7 a ~ p and {'ya-p, % :P '} = 
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2 8ab5 pp ' ■ The quantum Hall edge Hamiltonian now becomes 


^2? ^ hv J2 p H r b ~ V-pV-p) 


p> o 



p> o 


( 64 ) 



p> 0 



p>0 


Thus 



(65) 


p> 0 


up to a constant shift of the energy. This Hamiltonian is exactly two copies of a chiral 
Majorana Hamiltonian. The edge/domain-wall fermion Hamiltonian of the chiral p-wave 
superconductor will be 



( 66 ) 


p>0 


Finding gapless states on a domain wall of ji is an indicator that the phases with /x > 0 
and n < 0 are distinct. If they were the same phase of matter we should be able to 
adiabatically connect these states continuously. However, we have shown a specific case of 
the more general result that any interface between a region with /x > 0 and a region with 
/x < 0 will have gapless states that generate a discontinuity in the interpolation between 
the two regions. The question remaining is: Is /x > 0 or /x < 0 non-trivial? The answer is 
that we have a trivial superconductor for /x < 0 (adiabatically continued to /x —» — oo) and 
a topological superconductor for /x > 0. Remember that for now we are only considering /x 
in the neighborhood of 0 and using the continuum model expanded around ( p x ,p y ) = (0, 0). 
We will now define a bulk topological invariant for 2D superconductors that can distinguish 
the trivial superconductor state from the chiral topological superconductor state. For the 
spinless Bogoliubov-deGennes Hamiltonian, which is of the form 



(67a) 


p 


= (—2\A\p y , — 2\A\p Xl p 2 /2m — /x) , 


(67b) 


the topological invariant is the spectral Chern number defined in Eq. (11), which simplifies, 
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for this Hamiltonian, to the winding number 


rd) _ _ 

87 r 


d 2 p d ■ (d Pt d x d p d) = [ d 2 p d ■ (d Pi d x d p d) . (68) 


d I 3 


We defined the unit vector d = d/\d\, which is possible since \d\ ^ 0 due to the existence 
of a gap. This integral has a special form and is equal to the degree of the mapping from 
momentum space onto the 2 -sphere S 2 given by d\ + d\ + d\ — 1 . As it stands, the degree 
of the mapping d : M 2 —» S 2 is not well-defined because the domain is not compact, i.e., 
( PxiPy ) is only restricted to lie in the Euclidean plane (M 2 ). However, for our choice of the 
map d we can define the winding number by choosing an equivalent, but compact, domain. 
To understand the necessary choice of domain we can simply look at the explicit form of 

dip) 

^(p) = (~ 2 I A I Pm ~2|A|^,p 2 /2m - p) ^ 


\J 4| A| 2 p 2 + {p 2 /2 m — pY 

We see that limip^oo d(p) = (0, 0,1) and it does not depend on the direction in which we 
take the limit in the 2D plane. Because of the uniqueness of this limit we are free to perform 
the one-point compactiftcation of M 2 which amounts to including the point at infinity in our 
domain. The topology of M 2 U { 00 } is the same as S 2 and thus we can consider the degree 
of our map from the compactified momentum space (S' 2 ) to the unit d- vector space (S 2 ). 
Using the explicit form of the d-vector for this model, we find 

1 Arfe+x 


c' 1 * = ± 

7T 


d Z p~ 


2 p 2 +{£~h)' 


3/2 ' 


(70) 


The evaluation of this integral can be easily carried out numerically. The result is C^ 1J = 0 
for p < 0 and C (1 ) = 1 for p, > 0 , i.e., there are two different phases separated by a 
quantum critical point at /i = 0. Thus we have identified the phase which is in the chiral 
superconductor state to be p > 0 . 


1 . Argument for the existence of Majorana bound states on vortices 

A simple but rigorous argument can show us the presence of zero energy bound states in 
the core of vortices in a superconductor. Assume we have a chiral (p + ip) superconductor in 
two geometries: a disk with an edge and a cylinder with two edges. Since it is a topological 
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FIG. 7. Spectra of a Chiral superconductor in different geometries: (a) disk, (b) cylinder, and (c) 
disk with flux defect. Shown are the spectra of the chiral topological boundary modes including 
their finite-size quantization with level spacing v2n/L. If a n flux is inserted in the disk geometry 
(c), it binds an isolated zero-energy state. At the same time, a single zero-energy state appears on 
the edge. 


superconductor, the system will have chiral dispersing (Majorana) gapless inodes along the 
edges. In Fig. 7, the spectra are plotted versus the momentum along the edge, and they are 
qualitatively very different in the two cases. For an edge of length L , the smallest difference 
between two momenta along the edge is 2^t/L. The energy difference between two levels is 
v2i t/L, where v is the velocity of the edge mode. 

In a single particle superconducting Hamiltonian the number of total single-particle eigen¬ 
values is always even. This is clear from the fact that whatever the spinor of the nonsuper¬ 
conducting Hamiltonian is, when superconductivity is added, we have a doubled spectrum, 
so that every energy state at E > 0 comes with a counterpart at energy —E. When labeled 
by momentum quantum number, for a system with just one edge, like the disk, there cannot 
be a single state at momentum p — 0 at energy E — 0. If such a state was there, the 
spectrum would contain an odd number of states. Hence the spectrum of the linearized 
edge mode cannot have a state at E = 0, p = 0 on the disk. The one way to introduce 
such a state is to have antiperiodic boundary conditions, with the spectrum of the edge 
being at momenta 7 t(2 n + 1)/L, n G Z. On the cylinder, as two edges are present, periodic 
boundary conditions are allowed (as are antiperiodic, which can be obtained by threading a 
flux through the cylinder). 

We now add a single vortex inside the disk, far away from the edge of the disk. What is 
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the influence of the vortex on the edge? The vortex induces a phase 2i r in the units of the 
superconducting quantum he/ 2e, which means that the phase of A changes by 2n, and that 
of the electronic operators by 7r upon a full rotation around the edge. This implies that the 
antiperiodic boundary conditions on the edge without vortex changes to periodic boundary 
conditions in the presence of the vortex. The spectrum on the edge then is translated by 7 t/L 
compared to the case without the vortex, making it have an energy level at p — 0, E — 0. 
This would mean that the spectrum has an odd number of levels. However, this cannot be 
true, as we explained above, since the number of levels is always even. We are hence missing 
one unpaired level. Where is it? Since the only difference from the case with no vortex is 
the vortex itself, we draw the conclusion that the missing level is associated with the vortex, 
and is a bound state on the vortex. We also draw the conclusion that, since it is unpaired 
and really bound to the vortex, it has to rest exactly at E = 0, thereby showing that chiral 
superconductors have Majorana zero modes in their vortex core. 


2. Bound states on vortices in two-dimensional chiral p-wave superconductors 

Let us explicitly show that a vortex in a chiral superconductor will contain a zero 
mode. 12,13,15 For this calculation, which is a variant of our calculation for the existence 
of a Majorana mode at the interface between a topological and a trivial superconductor, we 
follow the discussion in Ref. 11. For this construction consider a disk of radius R which has 
H > 0 surrounded by a region with fi < 0 for r > R. We know from our previous discussion 
that there will be a single branch of chiral Majorana states localized near r — R, but no 
exact zero mode. If we take the limit R —>• 0 this represents a vortex and all the low-energy 
modes on the interface will be pushed to higher energies. If we put a 7r-flux inside the trivial 
region it will change the boundary conditions such that even in the R —» 0 limit there will 
be a zero-mode in the spectrum localized on the vortex. 

Now let us take the Bogoliubov-deGennes Hamiltonian in the Dirac limit (m —» oo) and 
solve the Bogoliubov-deGennes equations in the presence of a vortex located at r = 0 in 
the disk geometry in polar coordinates. Let A (r,8) = |A(r)|e ia(r b The profile |A(r)| for 
a vortex will depend on the details of the model, but must vanish inside the vortex core 
region, e.g., for an infinitely thin core we just need | A(0)| = 0. We take the phase a(r) to 
be equal to the polar angle at r. 
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The first step in the solution of the bound state for this vortex profile is to gauge transform 
the phase of A (r,9) into the fermion operators via T(r) —> e ialyr ^ 2 '^>{r). This has two 
effects: (i) it simplifies the solution of the Bogoliubov-deGennes differential equations and (ii) 
converts the boundary conditions of T(r) from periodic to anti-periodic around the vortex 
position r = 0. In polar coordinates the remaining single-particle Bogoliubov-deGennes 
Hamiltonian is simply 


BdG — ^ 

We want to solve 'Hedc'!' 


f 2|A(r)|e“(£ + i&) 

= E^> = 0 which we can do with the ansatz 


M) 


i 





MG 

|A(r')| 





( 71 ) 


(72) 


where J\f is a normalization constant. The function g(r ) is localized at the location of the 
vortex. We see that To(r, 9 + 27t) = — To(r, 9) as required. From an explicit check one 
can see that 'HBdG^o( r ; $) = 0. The field operator which annihilates fermion quanta in this 
localized state is 


7= j rdrdd ig(r) [—e l6 ^ 2 c(r, 6) + e l9 ^ 2 c\r,d)] , (73) 

from which we can immediately see that 7 = 7 !. Thus the vortex traps a single Majorana 
bound state at zero-energy. 


3. Non-Abelian statistics of vortices in chiral p-wave superconductors 

We have shown in the last Section that on each vortex in a spinless chiral superconductor 
there exists a single Majorana bound state. If we have a collection of 2 N vortices which are 
well-separated from each other, a low-energy subspace is generated which in the thermody¬ 
namic limit leads to a ground state degeneracy of 2 Ar . 16,1 ' For example, two vortices give a 
degeneracy of 2, which can be understood by combining the two localized Majorana bound 
states into a single complex fermion state which can be occupied or un-occupied, akin to the 
end states of the superconducting wire. From 2 N vortices one can form N complex fermion 
states giving a degeneracy of 2 N , which can be broken up into the subspace of 2 N ~ l states 
with even fermion parity and the 2 JV_1 states with odd fermion parity. As an aside, since 
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FIG. 8 . Illustration of the exchange of two vortices in a chiral p -wave superconductor. The dotted 
lines represent branch cuts across which the phase of the superconducting order parameter jumps 
by 27r. 


we have operators that mutually anti-commute and square to +1 we can define a Clifford 
algebra operator structure using the set of 2N 7 j. 

To illustrate the statistical properties of the vortices under exchanges we closely follow 
the work of Ivanov 14 and the discussion in Ref. 11. Let us begin with a single pair of vortices 
which have localized Majorana operators 71 and 72 respectively and are assumed to be well 
separated. We imagine that we adiabatically move the vortices in order to exchange the two 
Majorana fermions. If we move them slow enough then the only outcome of exchanging the 
vortices is a unitary operator acting on the two degenerate states which make up the ground 
state subspace. If we exchange the two vortices then we have 71 —> 72 and 72 —» 71 . However 
if we look at Fig. 8 we immediately see there is a complication. In this figure we have 
illustrated the exchange of two vortices and the dotted lines represent branch cuts across 
which the phase of the superconductor order parameter jumps by 27T. Since our solution of 
the Majorana bound states used the gauge transformed fermion operators we see that the 
bound state on the red vortex, which passes through the branch cut of the blue vortex, picks 
up an additional minus sign upon exchange. Thus the exchange of two vortices is effected 
by 

7i -» 72, 72 -> -71 • (74) 

In general, if we have 2N vortices, we can think of the different exchange operators Tij{j a ) 
which for our choice of conventions send 7 i 7 j, 7 j —)• — 7 j, and 7 —> 7 ^ for all k 7 ^ i, j. 
We can construct a representation of this exchange process on the Hilbert space by finding 
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a r(Tij) such that T{T l3 Yf a T 1 (T i j ) = T i3 { 7 a ). Such a representation is given by 

r(T ij ) = exp ^ (1 + 7j7i) • ( 75 ) 

Let us prove this by showing an explicit example for T 12 , which will have the transformation 
given in Eq. (74), 


(Tj 2 )7iT -1 (Tj 2 ) 

= \ (7i - 7 i727i + 72 - 7i) = 72, 

(76a) 

(Ti 2 )72T 1 (Ti 2 ) 

= ^ (72 - 7i + 727 i 72 - 72 ) = - 71 , 

(76b) 

(T 12 )73r- 1 (T 12 ) 

= 73'r(Ti2)r“ 1 (Ti 2 ) = 73 . 

(76c) 


Now that we have this representation we can illustrate the non-Abelian statistics. We 
start with four vortices with Majorana operators 71 , 72,73, 74 - To illustrate the action of the 
exchange operators on the four-fold degenerate ground state space we need to pair these 
Majorana operators into complex fermions 

a = ^(71 + 172), 

b = ^(73 + 174), 

The basis vectors of the ground state subspace can now be written as 

{ 10) a ® |0> 6 , | l)a ® |1) 6 , |l)a ® |0} 6 , |0) a 0 |1} 6 }, (78) 


a T = 2^ _ ^ 2 ), 
= ^(73 - m)- 


(77) 


where we have ordered the basis so that states of the same fermion parity are together. The 
notation |n) 0j & means, a)a \n) a = n\n) a and tfb \n)b = n|n)&. The set of statistical exchanges 
is generated by Tj 2 , T 23 , T 34 and we want to understand how these exchanges act on the 
ground state subspace. We can rewrite these three operators as 


t{Ti 2 ) = -^=(1 + 7271 ) = [! - i( aflt - « ta 

r(T 23 ) = —= [l — i (ba — ba) + tfa — &W)] 

V2 

t ( t m ) = -L [1 _ i(i*t _ 6 t 6) ] 

Taking matrix elements in our chosen ground state basis Eq. (78) we find 

/ (1 — i) 0 0 0 \ 


= ys 


0 (1 + i) 0 0 

0 0 (1 + i) 0 

\ 0 0 0 (1 - i) / 


(79a) 

(79b) 

(79c) 


(80a) 
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1 -i 0 0 
i 1 0 0 
0 0 1 -i 
0 0 i 1 

(1 + i) 0 0 0 

0 (1 - i) 0 0 

0 0 (1 — i) 0 

0 0 0 (1 + i) 

We see that with our basis choice T 12 and T 34 are Abelian phases acting on each state, while 
T 23 exhibits non-trivial mixing terms between the states with the same fermion parity. Thus, 
the form of T 23 represents non-Abelian statistics. Given an initial state |4in) = |0) a (8) |0)b, 
if we take vortex 2 around vortex 3 the final state is |-0f) = 4^ (|0) a 8) |0)b + i|l) a 8 |1)&). In 
principle one must also keep track of the Berry phase contribution to the statistical phase. 
Here we have only considered the wave function monodromy, however it can be proven 
that the Berry phase does not contribute in this case. The held of topological quantum 
computation is built on the idea that such exchange or braiding operations will lead to non¬ 
trivial quantum evolutions of the ground state which can be used for quantum computations. 

4- The 16-fold way 

We have now noticed that there are two characterizations of a topological superconductor, 
but they are seemingly different. First, the spectral Chern number is an integer C ( b £ Z. 
Directly related to it is the number of chiral Majorana modes on the edge, which in turn 
is related to an experimental observable, the thermal conductivity on the edge. Hence the 
system has a Z index, which becomes obvious when an edge exists. We then saw that a 
(p + ip) superconductor (i.e., a topological superconductor with Chern number equal to one) 
with a vortex threaded through it exhibits a Majorana zero energy mode at the core of 
the vortex. A (d + id) superconductor, with Chern number equal to 2 , would exhibit two 
Majorana modes in the core of the vortex. However, those two Majorana modes would be 
unstable towards single particle hybridization terms, which would push them away from 
zero energy, and leave the core of the vortex with no states in it. The generalization tells 
us that an even Chern number topological superconductor has no Majorana zero modes in 
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the vortex while an odd Chern number topological superconductor has one Majorana zero 
mode in its core. This shows that the defects (vortices) in a topological superconductor are 
classified by a Z 2 number (C^ mod 2). 

We now show that there is a third classification related to the idea of topological order . 18 
In the absence of an edge and in the absence of vortex defects, there is a Z 16 classification of 
topological superconductors indexed by C (1) mod 16, which can be put on solid grounds by 
the formalism of topological quantum held theory (TQFT) that we will introduce in the next 
Section. This shows that the edge-bulk correspondence needs revisiting - the bulk does not 
know if we add 16 edge modes or not, and hence that the edge contains more information 
than the bulk . 19 We first give a simple argument for the existence of a Zi 6 classification. 

We ask how we can classify the system in the absence of an edge. One way would be to 
compute the phases that wavefunctions can ackquire upon taking particles or quasiparticles 
around each other, ffowever, the system is made out of electrons (its a superconductor), 
so usually nothing special can happen to phases of electrons. The only “special” excitation 
of the superconductor is a vortex, so we will look at the phase that two vortices acquire 
upon exchange. We can calculate this with an argument. Take two copies of the (p + ip) 
superconductor governed by the Hamiltonian 

H = \ A jk{h,jh,k + 72j72,fc), (81) 

written in terms of Majorana operators 71 j for one copy and 72 j for the other copy. These 
operators can be combined into an complex fermion c 3 = (71 j + i 7 2 j )/2 in terms of which 
the Hamiltonian becomes 

H = \^A jk c]c k . (82) 

j,k 

This Hamiltonian has a “fake” U(l) symmetry given by our choice of Ajk for both Hamil¬ 
tonians. (Since the system is gapped, we expect our universal conclusions to hold even 
when this symmetry is stripped away). Thus, the system is a quantum Hall state of Hall 
conductance C (1) (in units of e 2 /h) if each of the superconductors had Chern number C (1) . 
We now ask what happens when we thread a superconducting vortex h/2e, which is equal 
to 7 r. Threading a flux 27r in a quantum Hall state of Chern number C (1) pulls C ( d electron 
charges to the vortex core through the Hall effect, hence a 7r flux pulls O 1 '/2 electron charges 
towards the core. We then try to compute the phase acquired when a vortex is exchanged 
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with another vortex. This is an exchange process, which is half a braid. A braid of two vor¬ 
tices is equivalent to C (1 )/2 electrons braided with a 7r vortex, giving rise to a phase 7rC (1 )/2 
upon a braid, and 7rC^/4 under exchange. Since this is the phase for exchange of vortices in 
two exactly identical superimposed superconductors, the phase for exchange in one of them 
is half that, 7rOb/8 = 27rC (1 yi6. This shows that the phase for vortex exchange is defined 
only mod 16. 

We will show this more rigorously within the framework of TQFT that we will introduce 
axiomatically in the next Section. Before doing so, let us summarize what we have learned 
about the vortices in chiral superconductors with odd Chern number in a language that 
anticipates the formalism that we will introduce. We have seen that well-separated vortices 
hold a Majorana zero mode at their core. When these vortices come together, the two 
Majorana modes hybridize and split, giving rise to two states which differ by their fermion 
parity. Let us call the Bogoliubov-deGennes vacuum 1 and the Bogoliubov quasiparticlc ip, 
and the Majorana fermion of the vortex cr. We can then formalize the fusion of two vortices 
by writing down a fusion rule 

cr x a = 1 + ip, (83) 

which basically tells us that combining two Majoranas can either go to a state with no 
fermion or at one with a fermion - the fermion parity (and density) would be different for 
the two states. Which one it is depends on the microscopies of the model. Hence a quantum 
state of two Majoranas has to be described by another quantum number, which describes the 
“fusion channel” of those two Majoranas - either the vacuum or the Bogoliubov quasiparticle. 
The fusion rule (83) allows for multiple fusion channels unlike the fusion rules (49) that we 
deduced for the toric code. This difference is a manifestation of the fact that the Majoranas 
are non-Abclian anyons, while the toric code anyons are Abelian. When two Bogoliubov 
quasiparticles fuse, they condense (form a Cooper pair) and go to the vacuum 

pj x ip = 1, (84) 

while the fusion of a Bogoliubov and a Majorana quasiparticle basically creates another 
Majorana 

ip x a = a. (85) 

This can be rationalized by thinking of the complex Bogoliubov quasiparticlc as made out 
of two Majoranas which then couple to the third Majorana. The Hamiltonian is a 3 x 3 
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antisymmetric matrix that necessarily has a zero eigenvalue which is another Majorana 
fermion coming as a result of the fusion. 


III. CATEGORY THEORY 

So far, we tried to gain some intuition about topologically ordered phases of matter by 
ways of examples. In this Section, we are going to define a framework that describes topo¬ 
logical order in 2D space in a unified and axiomatic way. At the same time, our description 
strips all nonuniversal details off the problem. A held theory with these properties is known 
as a topological held theory. It does not contain any information about energy scales of the 
problem. 

The topological held theory that we study is based on the mathematical concepts of 
category theory. 18,20,21 We will, however, try to keep the description as light as possible. For 
our purpose, we can view category theory as a generalization of group theory which is based 
on the fusion rules between anyon species that we have already encountered in examples. 
Consistent implementation of fusion defines a fusion category. Subsequently, we can impose 
more structure on the fusion category which elevates it to a braiding category, or a braided 
tensor fusion category. Our presentation will follow Refs. 18 and 21, while giving more 
examples of the use of the theory. 

A. Fusion Category 

A fusion category is based on a finite number of topological sectors (also called anyons, 
topological charges, or simply particles) which we will label 


( 86 ) 


a, b, c, • • • . 


For every charge a there exists a unique conjugate charge or anti-particle, that we denote 
by a. It is possible that an anyon is its own antiparticle a = a (for example the Majorana 


cr). There exists a unique vacuum sector denoted 1 (or sometimes 0). The fusion category 
is defined by its fusion rules 



C 
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where N c ah 6 Z + are nonnegative integers. We have already encountered two examples of 
fusion rules in the previous Section, namely the toric code with charges 1, e, m, and /, 


1 x e = e, 

1 x m — m, 

1 x / = /, 


e x m — /, 

e x / = m, 

m x / = e, 

(88) 

e x e = 1, 

m x m — 1, 

fxf = 1, 



and the so-called Ising anyon theory that we found realized by Majorana fermions in a chiral 
p -wave superconductor with charges 1, a, and T, 


1 x a — a, 1 x ip — ij) 

crxcr = l+^, a X -0 = cr, ^ X ^ = 1. 


(89) 


A principal difference between Eq. (88) and Eq. (89) is that the former has always only 
one fusion product on the righthand side, while the fusion of two a in the latter produces 
two outcomes. Hence, the x-product in the toric code can still be thought of as a group 
operation, while this is not possible in the Ising theory. We will see that this distinction 
coincides with the notion of an Abelian theory (toric code) and a non-Abelian theory (Ising). 

Does any choice of fusion rules, i.e., N% b G Z + , define a permissible fusion category? The 
answer to this question is negative, as we have to impose the following conditions on a fusion 
category. 


• The fusion rules be commutative 


a xb=bxa =* JVJ = JV£. 


(90) 


• The fusion rules be associative 

(a x b) x c — a x (b x c) => ]T = £ N L N S- (91) 

m m 

If we define the matrix N a with matrix elements ( N a ) bc = N% b , this relation becomes 
a vanishing commutator 

[Na, N c ] = 0, (92) 

which implies that all fusion matrices N a are diagonalized by the same eigenvectors. 
We will exploit this fact later. 
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• Fusion with the identity leaves any anyon unchanged 


a x 1 = a. 


( 93 ) 


The fusion product of a with its antiparticle a contains the vacuum with prefactor 1, 
i.e., 

«xa = l + ^ K- a c. (94) 

c^l 

There exists a solution the a consistency condition called pentagon equation , which we 
will discuss below. 


1. Diagrammatics 


Before we turn to the pentagon equation, we want to introduce a diagrammatic language 
that will facilitate computations within the fusion and braiding categories. In this formalism, 
we denote anyon a traveling forward in time as an upward oriented line. It is the same as 
the associated anti-particle traveling backward in time 



(95) 


A diagram with open anyon worldlines at the top and bottom represents a state in a 
Hilbert space that depends on the number and types of open anyon worldlines. A diagram 
without open worldlines represents an amplitude or complex number. The simplest nontriv¬ 
ial Hilbert space is the fusion space Vf b . Its dimension is given by the number of ways that 
anyons a and b can fuse into c, i.e., dim Vf h = Nf b . A basis in Vf b is denoted by 


(a, 6; c, p\ =: 


d a d b 


1/4 


a 


A 


e V: 


ab’ 


(96) 


b 


Here, /r = 1, • • • , Nf b labels the fusion multiplicity and the real positive prefactor [d c /d a d&) 1//4 
should be understood as a normalization constant at this stage. We adopt the normalization 
of Ref. 21. Likewise, we define the splitting space Vff b of the same dimension dim Vf b = Nf h 
and write a basis as 

dr 


I a,b\c,fi) =: 


d a d b 


1/4 “V 


b 


e V: 


ab 


(97) 
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Two propagating particles a and b live in a vector space V£ b = 0 C V^ b ®V^ b . The identity 


element I a b in V£ b is then represented by the completeness relation 


lab ^ ^ |^? /^) (^5 ^5 ^5 /^| 


(98) 


c,/x 


which we represent pictorially as 


a 


a 


E 


cL » 


V 


d a dbi .. 
c,H t l 


Ac 


(99) 


a 


A 


b 


The basis vectors in V^ b and VT furthermore satisfy the orthogonality relation 


(o, b, c, /r|a, 6, c , /i ) 


( 100 ) 


which we represent pictorially as 



In particular, we can choose c = 1 (a dashed worldline) to obtain 


( 101 ) 


a 




( 102 ) 


We can use this relation to determine the normalization constants d a , which are called 
quantum dimensions. We have the freedom to choose d\ = 1 and note that for all examples 
discussed here d a = d^. It follows that 



2. F-moves and the pentagon equation 

As noted above, we need to impose a further consistency condition to complete the 
definition of a fusion category. For this, we generalize the notion of associativity which we 
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imposed on the fusion coefficients by imposing associativity on the basis. The consistency 
equations basically say that observables only depend on the state of the particles at the 
beginning (the fusion channel) and at the end. Nothing in between can matter, up to phases 
and rotations in possibly degenerate spaces. We consider the Hilbert space of particle d 
splitting into three (not two) particles a, b, c 


Vt = E V? 0 V? 

e 


E v ? ® V T- 

5 


(104) 


There is hence a unitary transformation Ff ,c (“F-move”) between the two vector spaces 


a, 6; e, a) ® |e, c; d, j3) 


Y i F d bC ] k C ; d) ® k /; d, v) 

f,W 


(105) 


which reads diagrammatically 



(106) 


For a fusion category to be unitary , we require that (F^ 6c )f = (FJ ,,C ) 1 (as a matrix). Notice 
that F% bc is trivial if any of a,b,c = 1. 

The pentagon equation is the consistency condition diagrammatically represented in 
Fig. 9. It shows that there are two ways to build the same mapping between two vec¬ 
tor spaces out of F-moves. Each distinct solution, up to gauge freedom, of the pentagon 
equation 

F/% [K u ], k = E FFE [ff1„ ( 107 > 

h 

is a distinct fusion category (with the same fusion rules). Here and below we have suppressed 
the Greek indices that correspond to the fusion multiplicities. We will focus on theories of 
multiplicity 1 only, i.e., all N° b fusion coefficients will be either 0 or 1. As the F-moves relate 
different basis states, not all of them are gauge invariant (see below). However, there are 
some F symbols which are gauge invariant. Those are related to invariants of the theory 
called Frobenius Schur indicators. 
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a b c d 



FIG. 9. The pentagon equation defines a consistency relation that has to be imposed on a fusion 
category. 


3. Gauge freedom and its fixing 


Consider a gauge transformation on the basis states |a, 6; c) 


|a, b ; c)' = u ab |a, b ; c) 


(108) 


We only consider the case without multiplicities, i.e., Nf b = 0,1, for which uf b G C are scalars 
with \uf b \ = 1. Likewise, if N bc = 0,1 the F symbols are scalars with ^ = 1. In view 

of the definition Eq. (105), the F symbols are not invariant under the gauge transformation. 
They transform as 

Ff'LGGL (ms) 


1 e f n. a S q.bc 
U d Uj 


Furthermore, we need to set 


= Sic 


( 110 ) 


since the fusion of the identity particle can be added at any point in time to the worldline 
of any particle b without changing the state. 

From this, we can conclude that the following F symbols are gauge invariant 


m„= 


r pale] _ f pale] 

\ef~Y b d J ac > 




(Ill) 
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In a theory with no multiplicity, all these F symbols are equal to unity because they represent 
identity maps from spaces | b, c; d) into the same space 


= [Cii M = [FS'% = [fTL = [f?% = [ffl. = 1- (112) 


4- Quantum dimensions and Frobenius Schur indicators 


Having completed the definition of a fusion category, we now explore its structure. First, 
we shall properly define the quantum dimension d a of an anyon a that has already entered 
several relations as a normalization factor. Physically the definition of d a can be obtained 
from imposing isotopy invariance, which means the ability to remove bends in particle world¬ 
lines. This should be possible as long as lines are not crossed and end points are not moved. 
Bending a line slightly (so that the line always flows upward) is a trivial allowed move, but 
a complication arises when a line is bent so much that is acquires a turning point. The 
F-rnove associated with this type of bending is 



(113) 


Notice that the symbol [F““ a ] n is gauge invariant. Hence its value is a topological invariant. 
Since we know that the line in the left diagram is isotopically equivalent to a line going up, 
it should be, up to a phase, equal to a line going up. We conclude that [F“ aa ] n = Xa/d a , 
where Xa is a phase called the Frobenius Schur indicator. If a is its own antiparticle, Xa 
has to equal either +1 or —1. Since it can take different values, it is a topological invariant 
characterizing the fusion category. Interestingly, there are other Frobenius Schur indicators 
that characterize the theory as topological invariants. For example, one which we will not 
further elaborate on is connected to the trivalent vertex. 21 

We now know how to compute the quantum dimension through the F-syrnbols. However, 
there is another, easier way of computing the quantum dimensions. This is obvious once the 
space of states has been endowed with a completeness and an orthonormality relation, that 
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we gave in Eq. (99) and Eq. (101), respectively. One can use them to show the identity 


dadb ^ ^ Nab^ci 
c 


(114) 


for 



(115) 


Equation (114) is key to understand how the quantum dimensions follow from the fusion 
rules. It is again useful to render the fusion coefficients in a matrix form ( N a )b c = N£ b 
(b and c are the indices of the matrix N a ). Then, Eq. (114) is nothing but an eigenvalue 
equation for N a . We see that d a is an eigenvalue of N a and its eigenvector is the vector 
that contains all quantum dimensions d c . The existence of the real positive eigenvalue d a 
is a highly nontrivial fact. The Perron-Frobenius theorem, proved by Oskar Perron (1907) 
and Georg Frobenius (1912), asserts, in its weak version, that a real square matrix with 
nonnegative entries has a largest positive eigenvalue and that the corresponding (possibly 
degenerate) eigenvector has nonnegative components. We would like to use it to show that 
the eigenvalue d a is the largest eigenvalue of the matrix N a . By assumption, the vector with 
entries d c has only strictly positive components. Suppose we have another eigenvector v of 
N a with nonnegative components v c > 0 and eigenvalue /r a . Then the strict equality 


^2,Vcdc> 0 (116) 


holds since all d c are strictly positive and at least one v c is strictly positive as well. From 



= 22 d bKb V b = Va 

b,c 


y 4 ,'!-,. 


(117) 


it thus follows that the eigenvalues d a and fi a are equal. In other words, any nonnegative 
eigenvector of N a has the same eigenvalue d a . This includes the eigenvector of the largest 
eigenvalue of N a , which is nonnegative due to the Perron-Frobenius theorem. Hence, d a is 
the largest eigenvalue of N a . 
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Using the fact that d a is the largest eigenvalue, we can now give a more physical inter¬ 
pretation of the quantum dimension. For that, consider the fusion of some anyon a with 
itself n times 

Q X Q X ' ~ 'X CL y ^ (-^a)aci (-^a)ciC2 X • • • X (No) Cn-2Cn_i C n _i- (118) 

71, Cl ,C2 }••* ,Cn — 1 

The righthand side contains the (n — l)-th power of the matrix N a . Hence, approximating 
N a by its highest eigenvalue, we conclude the the dimension of the fusion space of n anyons 
of type a is dominated by the quantum dimension dim((J) c V£ a ) ~ d ™ for large n. In other 
words, the quantum dimension tells us how fast the Hilbert space for a particle grows! Any 
non-Abelian particle has a quantum dimension strictly larger than unity (Abelian particles 
have quantum dimension unity). 

5. Examples 

Before moving on to impose more structure on the fusion category in order to obtain a 
braiding category, we shall briefly follow up on the two examples of semion TQFT and Ising 
TQFT. 

Semion TQFT — The simplest nontrivial TQFT has one particle s besides the identity 
and the semion fusion rules 

s x s — 1, 1 x s — s. (119) 

The theory is Abelian, i.e., d s = 1. Let us solve the pentagon equation for this theory. 
There is only one F-symbol, which is not entirely determined by gauge fixing alone, namely 
[Fr ], j. We can deduce its allowed values from the pentagon equation 

[*?"]., Fr 1 ],, = [friii [FTLiFri,, ( 12 °) 

which, using Eq. (112), yields the two possibilities 

[FT],, = ±1. (121) 

In fact, [F, sss ] n is equal to the Frobenius Schur indicator mentioned above and the two 
values ±1 distinguish two different fusion categories. The choice +1 is trivial, while the — 1 
is what is commonly called the semion TQFT. For example, it is realized in the u — 1/2 
Laughlin state of bosons in the fractional quantum Hall effect. 
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Ising TQFT — We gave the fusion rules of the non-Abelian Ising TQFT in Eq. (89). We 
now want to use Eq. (114) to compute the quantum dimensions of the anyons. For that we 
note that the fusion matrix [in the basis (1, cr, /)] of the a anyon is given by 


N„ = 


'0 10 ' 
1 0 1 

V 010 / 


( 122 ) 


Its eigenvalues are given by ±a/2 and 0, the largest of which = \/2 is the quantum 
dimension of the a particle. Its corresponding eigenvector (d\, d a , d^) = (1, a/ 2, 1) indeed 
contains the quantum dimensions of all anyons. The quantum dimension y/2 is compatible 
with our explicit calculation of the degeneracy resulting from Majorana bound states in the 
vortices of a p-wave superconductor. For example, 2 vortices with one Majorana state each 
gave rise to a degeneracy a/2 = 2 of the state. 


B. Braiding Category 


A bare fusion category has no means to relate the two fusion spaces and V ba . The 
physical operation that corresponds to such a relation is an exchange between the particles 
a and b. In a braiding category, exchange is a map R a b ■ —> V ba . A double exchange 

is an automorphism in a given fusion space Rb a R a b '■ V/ 6 —> V/ 6 . For fusion processes 
without multiplicities, for which the space is one-dimensional, Rb a R a b is thus a phase (it is 
represented by a matrix if the multiplicity N° b is larger than 1). Upon braiding the state 
changes by flipping the particles, by convention, 


Rab | a,b;c,n) = ^ \b,a; c, u). 

V 


Diagrammatically, this operation (“/?, move”) is represented as 



(123) 


(124) 


in the case without multiplicities. 

For this braiding relation to define a consistent braiding category, the R and F symbols 
have to satisfy two consistency relations called hexagon equations. (The two hexagon equa¬ 
tions differ only in the directions of the braids involved; one uses R, the other one uses R~ 1 . 
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FIG. 10. The hexagon equation defines a consistency relation that has to be imposed on a braiding 
category. 


For simplicity, we only present the equation for R here.) As with the pentagon equation for 

the F symbols alone, the hexagon equation represents a way to mix F and R moves to get 

between the same two diagrams in two different ways. Analytically, it reads 

K b = [Ft\e f R°J [Ci “] /9 (125) 

and they are shown diagrammatically in Fig 10. The matrix R^ b is unitary, and if any of a 
or b is the identity particle then the R matrix is unity (braiding with the vacuum is trivial) 

Rl b = Rf = 1. (126) 


1. Topological spin 

Every quasiparticle type a in a braiding category carries another topological quantum 
number besides its quantum dimension, namely the complex-valued topological spin 6 a with 
0 a | = 1. Physically, the topological spin can be interpreted as the phase resulting from a 
rotation of the particle around its own axis by 2n. It is —1 for fermions, but can be fractional 
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for other particles. Diagrammatically, we can define 6 0 as 



(127) 


We bend the world line of a particle a, bring it down, then bring it back to itself, this 
time braided, then take it to infinity. This corresponds, np to a phase, to the configuration 
without bending. 

To see how the topological spin is related to the R symbols, we want to evaluate a 
diagram that represents an amplitude, i.e., it is closed. For that, we simply connect up the 
open worldlines in Eq. (127) with an a worldline and obtain the definition 



(128) 


We evaluate this diagram by inserting an identity and applying an R move (by the definition 
of the R symbols, exchange is only defined between worldlines that are in a definite fusion 
channel) 


a. 



a = 



E dc R' 

\!e R ‘ 


aa'f a 


a'' 



(129) 


Combining Eqs. (128) and (129), we obtain the definition of 9 a in terms of the R symbol 




r\ 


(130) 


Here, the matrix trace Tr is relevant for theories with multiplicities, in which case i?““ is a 
square matrix with dimension equal to the multiplicity of the fusion channel c of a x a. 
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Another identity, which requires bending to be proven, is 


R 


aa 

1 


0* a Xa 


(131) 


where \a is the Frobenius Schur indicator. 


2. Ribbon equation 


We are now in good shape to prove an important identity relating topological spin to the 
R matrices that is called the ribbon equation 


r>ab T)ba 
it n lt n — 


0 a 0b 


(132) 


where / is the identity element in the space V^ b . Physically, it equates the operation of 
twisting each worldline in a splitting diagram by 27T with the operation of braiding the 
splitting products. We can prove the ribbon equation via diagrammatic manipulations. For 
that, we evaluate the same diagram in two different ways. For one, we can use the topological 
spin of the split particle c 



(133) 
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On the other hand, we can use a combination of topological spins and R moves 



(134) 


Together, Eq. (133) and Eq. (134) yield the ribbon equation (132). 


3. Vafa’s Theorem 


An important theorem relating the topological spin and the structure constants N° b is 
Vafa’s theorem. 22 It shows that the topological spin is a rational number. We will not derive 
it here but point the reader to Ref. 18 for an easy derivation. Vafa’s theorem proceeds by 
writing down in matrix form the two hexagon equations, one for R and one for i?” 1 , dividing 
them and using R“ b R ba = 0 c /(0 a 0b)I to obtain 


n 


6 a 0 b 


n 

/ 


d f \ N bf N ad 


0 a 0d 


n 


Ondr 


(135) 


As an example, we will later use Vafa’s theorem to deduce the spin of the a particle in the 
Ising TQFT. 


C. Modular matrices 

In the last Sections, we have constructed the fusion and braiding category from the F 
and R moves and deduced the universal data (quantum dimensions d a and topological spins 
6a) that characterizes the anyons in this category. 

For TQFT in physical systems, i.e., quantum liquid ground states of matter, it is impor¬ 
tant to address how the universal information about the TQFT can in general be accessed. 
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We know that topological properties cannot be accessed by local measurements. In contrast, 
it turns out that there is a set of global measurements that can be used, namely the action of 
automorphisms on the manifold over which the system is defined. Automorphisms are trans¬ 
formations that map the manifold back to itself and form the so-called mapping class group 
of the manifold. Here, we will explore the most standard case, namely (2 + l)-dimensional 
systems with periodic boundary conditions, in which case the manifold is a torus. Automor¬ 
phisms on the torus form the modular group, which has two generators. The first generator 
S exchanges the two coordinate axes. The second generator T changes the angle between 
the coordinate axes. 

A TQFT on the torus exhibits a topological ground state degeneracy, where the number 
of ground states is equal to the number of anyons in the theory (including the identity). It is 
the representation of the S and T operations in this ground-state manifold that reveals infor¬ 
mation about the nature of the TQFT. We call these representations the S and T matrices. 
Instead of deriving the S and T matrices from the action of the respective transformations, 
we will simply give their definitions within the TQFT here and subsequently explore their 
properties (see Ref. 23 for a more complete discussion of the connection). 


1. The S matrix 


In the TQFT, the S matrix is defined diagrammatically as 



where D is the total quantum dimension of the TQFT 



(136) 


(137) 


(We choose the same convention for the diagrammatic definition of S a b here as in Ref . 11 
which is different from the one chosen in Ref. 18 ) We can relate it to the topological spins 
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and fusion coefficients via the following diagrammatic manipulations 




(138) 


From here, we use the ribbon Eq. (132) and Eq. (136) to arrive at the final expression for 
the S matrix 


Sa„ = 1 E N * $C 


0 a 0b 


dr. 


(139) 


2. Verlinde Formula 


We now derive a fundamental formula in both TQFT and in conformal field theory. This 
formula relates the 5-matrix to the fusion coefficients. It allows us to find the set of braiding 
phases among the anyons, but not always the topological spin of every particle. 

We start by showing that 

Si, = % (140) 

To prove this, we first observe that 



(141) 


holds, because we can close the x worldlines in this diagram to obtain amplitudes and an 
identity loop can be added without changing the value of the diagram. Relation (141) can 
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be used to obtain the line of equalities: 


DS ax = 




(142) 


which yields Eq. (140). Equation (140) says that the first column of the S matrix contains 
only positive numbers, the quantum dimensions of the theory. 

The next step is to derive the all-important Verlinde formula. Staring from two copies of 
Eq. (141), we can perform the following set of diagrammatic manipulations 



£ K„C> = 


^lx 


Ai X 


c 


X. 


yielding the Verlinde formula 

y,N^S CI = S bI §=. (144) 

Oi r 
c 

This is a remarkable equality, which takes again the form of an eigenvalue equation for the 
fusion matrices N a . We had already encountered an eigenvalue equation for N a when solving 
for the quantum dimensions, i.e., the Erst column of the S matrix. The Verlinde formula 
says that the S matrix contains both the eigenvectors and the eigenvalues of the ( N a )b c 
matrices. If the S matrix is unitary (and thus has an inverse), we call the theory a unitary 
modular category. In this case we have a simple relation between the fusion coefficients and 
the S matrix 


A'„ = SD.S-\ 


(145) 


where we have dehned the diagonal matrices 


\^a)mn 3 am /S\ m 5 rnn . 


(146) 
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This also implies that all the N a matrices are commuting [which we already knew from 
Eq. (92)], as they are diagonalized by the same eigenvectors. 


3. Obstruction for theories with multiplicities 


In this Subsection, we want to illustrate how we can use the structure of the braiding 
category to discard certain fusion categories as unphysical (or at least not unitary). As an 
example, we ask which theories are possible with only two particles, the identity 1 and a 
particle s. For s to have an inverse, all possible fusion rules 


sxs=l+ms 


(147) 


are labelled by a nonnegative integer m. The case m = 0 is the semion TQFT of Eq. (119). 
For m — 1 we have the non-Abelian Fibonacci fusion rules. 

We would like to answer the question for which m > 1 we can define a consistent modular 
braiding category. Observe that the fusion matrix 



yields the quantum dimension 


da 


m 


+ \/m 2 + 4 


From the definition of the S matrix, we have that 
S.. = jE Ki= hd, = S u] Sn = h 


—(1+mM, 


(148) 


(149) 


(150) 


It is already possible to see that something will go wrong for large enough m if we demand 
a unitary theory, i.e., a theory with a unitary S matrix. For a unitary matrix, all the matrix 
elements have to be less or equal to 1. When m is large, d s is proportional to m (and so 
is D) but the S ss matrix element has a md s /D which is proportional to m in the large m 
limit. Hence we forsake unitarity. 

Lets us calculate the exact m where unitarity breaks down. Imposing unitarity of the S 
matrix yields 


S ] S 


1 


D 2 


d s 



1 +mO*d s \ 

0*sf ) 


d s (1 + 


1 -\-m9 s d s 

0 ? 


d 2 + (1 + m6*d s )( 1 + mOgdg) 


(151) 
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and we find 


(1 + m6*d s ){ 1 + m6 s d s ) = 1, 


(152) 


which reduces to 


0 S + 9* S + md s = 0, 


(153) 


where we have used that the value of d s is positive. Now, md s is a positive number growing 
with m, while 6 S is a phase so the sum with its conjugate cannot be smaller than —2. Hence 


md s <2 —>■ m 2 + mVm 2 + 4 < 4 


(154) 


with the solutions m = 0 and m, — 1. We can see that already m = 2 gives a left side too 
large. We conclude that the semion and the Fibonacci TQFT are the only allowed modular 
unitary theories with one nontrivial anyon s. 

4■ The T matrix 

The T-matrix is diagonal and simply given by 


(155) 



D. Examples: The 16-fold way revisited 

To motivate our study of TQFTs and as an example of topologically ordered phases, we 
have studied in Sec. IIB 4 Kitaev’s 16-fold way of classifying topological superconductors as 
gauge theories from their bulk properties. Equipped with category theory understanding of 
TQFTs, we now want to revisit this classification, as it provides us with several examples 
of TQFTs. In particular, we want to characterize the theories from their S and T matrices. 

1. Case: Q 1 ) odd, 

If we couple an odd number of layers of spinless chiral p- wave superconductors, the core 
of each vortex still carries an unpaired Majorana state. For that reason, the vortices will 
have Ising fusion rules (89) in this case. Using the Verlinde formula, we can compute the S 
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matrix in the basis (1, a, 0) 


( 


1 


\ 


V 


S — y /\2 0 -a/2 • (156) 

1 -y/2 1 y 

If, in contrast, we compute the S matrix via its dehnition (139), we find in particular for 
the (cr,cr) matrix element 




l + ^i/ 

2 91 


(157) 


For this to vanish, we conclude that 9^ = —1, that is, / is a fermion. We cannot obtain 9 a 
from similar relations using only the S matrix. Rather, we can use Vafa’s theorem (135) to 
deduce the spin of the a particle. If we take a = b = d = e = a we find 


n< 


Oc TT / Of 


ej n 


w (7 W CT 


9 n 


9Jr 


or 


o± = JL 

01 ei 9191 9 / 


which gives 


Since 9^ = —1, we have 


01 = 0}. 


(158) 

(159) 

(160) 


el = e t = -l, (lei) 

i.e., the phase of the topological spin of the Majorana is an odd-integer multiple of 1/16. 
We can thus discriminate 8 different TQFTs with Ising fusion rules by the values 

cd) 


0„ = e 2 "-nr, 


(162) 


where C ( b is the [odd, as needed by Eq. (161)] Chern number or the number of stacked 
spinless chiral p-wave superconductors. 


2. Case Cf 1 ) = 2 mod 4 

If we stack C-k = 2 mod 4 layers of chiral p -wave superconductors, we had argued before, 
heuristically, that the system can be described as one species of spinlcss Dirac fermions 
with Chern number C (1) = C^/2. In this system the 27T flux i/j binds odd integer charge 
C (1 \ that is, the 27r flux is a fermion -0. However, unlike in the quantum Hall effect, the 








superconducting n or (—vr) fluxes are allowed topological excitations that bind half-integer 
charge. Let us denote the ir flux with charge 0^/2 by a. Fusing such a it flux with the 
fermion (or 2tt flux) gives another excitation, a (— 7 r) flux, with 3C ( h/2 charge that we call 
a. (We note that 47 t flux is identified with zero flux, as this corresponds to a charge 2C (1 ^ 
object, i.e., C (1 ) Cooper pairs that can be absorbed by the condensate.) Two fluxes of either 
type a or a thus fuse into ijj. This motivates the following fusion rules 


a x a = 1, a x ?) = a, a x if) = a, 


axa = axo = 'i), ^ x ^ = 1. 


(163) 


These fusion rules are Abelian, so that the quantum dimensions are d\ = d a = da = = 1. 

Given the fusion matrices, we can compute the S matrix as the matrix of their simultaneous 
eigenvectors. In the basis it can take one of two forms 



fl 1 1 l\ 


(ll 1 1 \ 

S C) = i 

1 -i i -1 

. S' 2 ' = 1 

1 i -i -1 

2 

1 i -i -1 

2 

1 -i i -1 


1 

1 


1 

1 


From S aa = Saa = ±i/2, we conclude 


9 


2 

a 


Tb 


(164) 


(165) 


while S a ^ = Sai); = - 1/2 gives 


9 


a 


9 a, 


(166) 


with the help of Eq. (139). We conclude that there are four possible theories with topological 
spins 


9 a = 9~ a = e 


27ri 


C (!) 


(167) 


3. Case (T 1 ' = 0 mod 4 

If we couple 4 layers of chiral p -wave superconductors, the it superconducting vortices 
(let us denote them e) bind a full electron charge. However, the electrons exist also as free 
fermionic quasiparticles / in the theory. Hence, the fusion of e with / should yield a new 
excitation m that is a 7r vortex stripped of its charge. Together, e, m, and / obey the toric 
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code fusion rules Eq. (88). The S matrix of the theory in the basis (1, e, m, /) can take one 
of two forms 


S<-> = i 

fl 1 1 l\ 

1-11-1 

. s (21 = l 

fl 1 1 i\ 

11-1-1 

2 

h ^ 

h ^ 

1 

h ^ 

1 

h ^ 

2 

h ^ 

1 

h ^ 

h ^ 

1 

h-L 


1 

h 

1 

h ^ 

h-L 


h-L 

1 

h ^ 

1 

h ^ 

h ^ 


(168) 


If we assume that 9f = —1, we find from that 9\ = 9 1 2 3 4 5 6 * 8 9 10 11 12 m = 9 e 9 m = — 1, while the theory 
with S ( 2) has 9\ = 9 2 m = 9 e 9 m = 1. In total, we have four possibilities 

9 e = 9 m = e 2m (169) 


with C (1) = 0mod4 as allowed theories. This concludes Kitaev’s 16-fold way. 18 
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